Junctions of three quantum wires for spin 1/2 electrons 
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We study the effects of electron-electron interactions on the transport properties of a junction 
of three quantum wires enclosing a magnetic flux. The wires are modeled as single channel spin- 
1/2 Tomonaga-Luttinger liquids. The system exhibits a rich phase diagram as a function of the 
electronic interaction strength, which includes a chiral fixed point with an asymmetric current flow 
highly sensitive to the sign of the flux, and another fixed point where pair tunneling dominates, 
similarly to the case of spinless electrons. While in the case of spinless electrons the perturbations 
that correspond to unequal couplings between the three wires are always irrelevant, we find that, 
when the electron spin is included, there are small regions in the phase diagram where a current 
flows only between two of the wires and the third wire is decoupled. 
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I. INTRODUCTION 

The transport properties of quantum wire systems 
have been the subject of intensive investigation in the 
past one and a half decade, both because they can have 
practical applications in nano-electronic circuits and be- 
cause they provide an experimentally realizable way for 
understanding exotic properties of one-dimensional inter- 
acting electron systems. As opposed to two- and three- 
dimensional electronic systems, onc-dimcnsional (ID) 
systems cannot be described by Fermi Liquid theory; in- 
stead, they are described as Tomonaga-Luttinger liquids 
(TLL) 0, 0, IE El , and experiments on carbon nanotubes 
are consistent with this description [E [f| . 

Recently, there has been a number of studies of junc- 
tions of multiple TLL wires_0Jl, JE EE El El EE E3, 

EE EE EE EI EE HE HI, HaM HlHi ■ such studies 

are relevant because junctions of three or more quan- 
tum wires would inevitably appear in any quantum cir- 
cuit. New tools and methods for investigating junctions 
of three quantum wires with spinless electrons were pro- 
posed in Rcf. Q. These new methods allowed for the 
identification of a low energy chiral fixed point with an 
asymmetric current flow that is highly sensitive to the 
sign of the magnetic flux enclosed at the junction. There 
are, however, important outstanding issues in the three- 
wire junctions problem, for instance, the behavior of the 
more realistic model in which the electron spin is taken 
into account. 

Even in the case of tunneling between two wires, the 
inclusion of spin degrees of freedom already brings about 
a rich phase diagram in the charge and spin interaction 
parameter space [IE HE EH- For example, one can find 
a situation where the charge conductance vanishes while 
the spin conductance does not, or vice versa. In the case 
of junctions of three quantum wires (Y-junction) for elec- 
trons with spin, the phase diagram becomes much richer. 
In addition to a phase similar to the two-wire case in 
which electron pair tunneling dominates while the spin 
conductance vanishes, or vice versa, the chiral fixed point 
revealed previously in the spinless Y-junction case per- 



sists when the spin degrees of freedom are taken into 
account. Moreover, while the strong asymmetric limit in 
which one of the wires is totally decoupled was proved al- 
ways unstable for the spinless case Q , we find a small re- 
gion in the coupling parameter space in which this asym- 
metric fixed point is stable. 

The paper is organized as follows. In Sec. [II] we sum- 
marize the results and present the phase diagram in cou- 
pling constant space at zero temperature. We also dis- 
cuss the conductance tensor corresponding to each stable 
fixed point. In Scc. lIIII we present our effective model for 
the junctions of three quantum wires for spin-1/2 elec- 
trons. In Sec. IIVI we review the Delayed Evaluation of 
Boundary Condition (DEBC) method introduced previ- 
ously for identifying the stable fixed points for the Y- 
junction system with spinless electrons. Taking advan- 
tage of charge/spin separation in one dimension, we gen- 
eralize this approach to the system with spin-1/2 elec- 
trons and illustrate the method by describing the phase 
diagram of a junction of two quantum wires. In Sec. [V] 
we apply the DEBC method described in Sec. IIVI to the 
Y-junction with spin 1/2 electrons and determine the sta- 
bility of each fixed point. In Sec. lVII we briefly review the 
method of Boundary Conformal Field Theory (BCFT) 
and apply it, as a warm up, to study a junction of two 
quantum wires. We then apply BCFT methods to the 
Y-junction problem in Sec. IVIII and obtain results con- 
sistent with those found with the DEBC method. We 
briefly conclude in Sec. IVIlfl 



II. SUMMARY OF RESULTS 

In this section, we summarize our study of the stability 
of the fixed points associated with the different boundary 
conditions of Y- junctions for spin-1/2 electrons and their 
corresponding conductance tensors. Based on the stabil- 
ity of each phase, we propose a zero-temperature phase 
diagram in terms of the TLL interaction parameters, g c 
and g s . The geometry of the device is depicted in Fig. [U 
where three identical quantum wires are attached to a 
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FIG. 1: [Color online] A junction of three quantum wires with 

a magnetic flux threading through the ring. The J 12 3 are the 
charge and spin currents arriving at the junction from each of 
the three wires. The lighter line (Blue online) indicates the 
charge sector, and the darker line (Red online) presents the 
spin sector. 



ring, with equal couplings and thus having a Z^ symme- 
try. The ring can be threaded by a magnetic flux and 
the quantum wires arc modelled as TLLs characterized 
by the interaction parameters g c and g s . We will only 
consider spin conserving transport. 

Due to the charge and spin separation in the bulk of 
TLL quantum wires, independent boundary conditions 
can be imposed on the charge and spin degrees of free- 
dom. Hence, we introduce the notation B C B S in which B c 
and B s represent the corresponding boundary condition 
in the charge and spin sector respectively. For instance, 
DN represents Dirichlct boundary condition (B.C.) in 
the charge sector and Neumann B.C. in the spin sec- 
tor. Moreover, since each combination of the boundary 
conditions corresponds to a fixed point in the framework 
of the renormalization group (RG), the fixed point cor- 
responding to B C B S BC will be referred to as B C B S fixed 
point. The physical meaning of each boundary condition 
will be discussed in this section and the detailed analysis 
and assumptions of the model are given in later sections. 

The charge and spin conductance tensors associated 
with each fixed point are important physical response 
functions characterizing the Y-j unction for spin 1/2 elec- 
trons. Within the linear response theory, the total cur- 
rent Ij flowing into the junction from wire j is related to 
the voltage applied at the wire k through 



T 3 = /\2 G jk V k> 



(2.1) 



where j, k = 1,2,3 and Gjk is the 3x3 conductance 
tensor. One can similarly define, 



(2.2) 



where A^ fe = fi\ — (if. is the chemical potential differ- 



ence between up and down spins in lead k, and e is the 
electron charge (with this definition, spin and charge con- 
ductances are both measured in units of e 2 /h). 



Note that current conservation implies that: 

0. (2.3) 
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Furthermore, a common voltage applied to all three wires 
results in zero current. Thus: 



(2.4) 



For Z$ symmetry junctions, the conductance tensor 
takes the form Q 
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(2.5) 



where we separate the symmetric and anti-symmetric 
components of the tensor, and G$ and Ga are scalar con- 
ductances. (The 6ij are defined as follows: £12 = £23 = 
£31 = 1, £21 = £32 = £13 = -1 and en = e 22 = £33 = 0.) 
The anti-symmetric component G A will only appear 
when time reversal invariance is broken. Even in the pres- 
ence of the magnetic flux which would normally break 
time reversal symmetry (TRS), Ga may vanish at some 
low energy fixed points, in which case time reversal sym- 
metry is restored. However, in the absence of the Z3 sym- 
metry (in the asymmetric fixed point), the condition of 



Eq. (|2.5[) becomes unnecessary. Observe that G s = Gjj 
represents the conductance of each wire when zero volt- 
age is applied to the other two wires (i.e., a potential 
difference is applied between one of the wires and the 
other two, which are held at the same potential). 

We are mostly interested in the stable RG fixed points 
which describe the physics in the low energy limit, i.e. 
low voltage bias and low temperature. As in the case of 
junctions of two quantum wires, the interaction parame- 
ters control the stability of the fixed points and determine 
the phase diagram. Below we list our results for con- 
ductance tensors and for the basins of attraction around 
each fixed point (with their corresponding boundary con- 
ditions) as function of the interaction parameters, g c and 
9s- 



A. Neumann BC in both charge and spin sector 

The NN boundary condition corresponds to a fixed 
point in which the hopping amplitudes between the wires 
are zero and the three quantum wires are totally decou- 
pled from each other, as illustrated in the inset of Fig. [5] 
When the interactions in the quantum wire are repul- 
sive, the hopping amplitudes decrease alon g th e RG flow 
both injunctions of two quantum wires [261.1 2 71 1 2 81] and in 
junctions of three quantum wires for spinless electrons 0]. 
Generically, the window of stability for a fixed point with 
decoupled wires should be independent of the number of 
wires in the junction. Notice that the effective flux <j> has 
no effect in the decoupled limit. 
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FIG. 2: [Color Online] The painted area (Red) shows the 
attractive basin of NN fixed point. The border shows the 
marginal line of scaling dimension A = 1 for all leading order 
perturbations. The inset depicts the physical consequence of 
the fixed point related to the NN BC, a disconnection for both 
charge and spin sectors. 



The three most relevant perturbations at the NN fixed 
point are: 

(1) Single electron hopping between two wires, with 
scaling dimension Ann = ^(— + — )■ 

(2) Electron pair singlet hopping between two wires, 
with scaling dimension Ann = —■ 

(3) Exchange of electrons with opposite spins between 
two of the wires (thus carrying a spin current) , with scal- 
ing dimension Ann = j-- 

By requiring that all scaling dimensions corresponding 
to these leading order perturbations Ann > 1, we can 
identify the attractive basin in the interaction parameter 
space, as shown in Fig. O The conductance tensor for 
the NN fixed point is zero for both the charge and spin 
sector. 
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FIG. 3: [Color Online] The painted area (green) shows the 
attractive basin for the DD fixed point. The border shows 
the marginal line when all leading order perturbations have 
A = 1. The inset shows one of the tunneling processes as- 
sociated with the DD fixed point. The conductances in both 
charge and spin sectors are enhanced by the pair hopping (or 
Andrcev-like) processes in both spin and charge sectors. 



presence of the effective flux (j> which breaks TRS has 
no physical consequence and the TRS is restored in DD 
fixed point. 

We shall study the leading order perturbations to de- 
termine the stability of the DD fixed point. Since several 
boundary operators possess the same scaling dimensions, 
we only list the dimensions of the leading order pertur- 
bations without specifying the corresponding operators, 



A 1 / , . n 

A DD = glflc + .9s), A DD = — 



A 3 - 

^DD — 



2g s 



(2.8) 



The attractive basin of the DD fixed point, shown in 
Fig. [3j can be obtained by requiring the dimensions of 
these leading order perturbations to be larger than one. 



ND and DN BC 



B. Dirichlet BC in both charge and spin sector 

The DD BC corresponds to a fixed point where 
the paired-electron hopping and Andrecv reflecting pro- 
cesses, depicted in the inset of the Fig. [3j dominate the 
tunneling processes in both charge and spin degrees of 
freedom. Hence, the charge and spin conductance will 
be enhanced. Moreover, the conductance tensor takes 
the symmetric form in which the scalar conductance 



G 



c(s) 



2 x (4<? c ( s )/3)(e 2 //i) and the antisymmetric com- 



ponent Ga — 0, and is given by 

2.9c( 



G% S) \ DD =2X^(%-1) 



(2.7) 



Here, the factor of 2 comes from the doubling of the de- 
grees of freedom due to the spin. Since Ga — 0, the 



We have seen that the charge and spin current termi- 
nate or flow at the same direction for both NN and DD 
fixed point. However, phases with decoupled charge and 
spin degrees of freedom are also possible. Indeed, the ND 
and DN BC correspond to fixed points where the charge 
or spin degrees of freedom are disentangled at the bound- 
ary. The ND fixed point possesses a pure spin current; 
likewise, the DN fixed point possesses a pure charge cur- 
rent. The dominant processes corresponding to these two 
fixed points are illustrated in the insets of Fig. |U Notice 
that similar phases also exist in the system of junction of 
two quantum wires for spin-1/2 electrons. 

Observe that the dominating process of the ND fixed 
points carries no net charge current while that of DN 
fixed point carries no net spin current, hence the charge 
and spin conductance vanish at the ND and DN fixed 
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FIG. 4: [Color Online] The painted area shows the stable 
region for the ND and DN fixed points. The borders (Blue 
for ND BC and Yellow for DN BC) show the marginal line 
of scaling dimensions A = 1 for all leading order perturba- 
tions. The inset shows a pictorial representation of the ND 
and DN fixed points. The fundamental process associated 
with the ND fixed point leads to an enhancement of the spin 
conductance due to a spin exchange between two wires, while 
the fundamental process associated with the DN fixed point 
leads to an enhancement of the charge conductance due to a 
electron-hole exchange. Notice that the processes described 
here occur between any two arbitrary wires. 



point, respectively. The conductance tensors, taking the 
symmetric form with Z% symmetry, are given by 

G%\nd = 2^ (36 jk - 1), G c jk \ ND = 0, (2.9) 

G%\dn = 2^ (35 jk - 1), G' Jk \ DN = 0(2.10) 

The leading order perturbations of the ND BC have 
the dimensions, discussed in detail in Sec. IVland IVII1 



1 
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The basin of attraction relative to the ND fixed point is 
obtained by requiring all A J ND > 1 and is depicted in 
Fig. [4] By exchanging g c and g s in Eq. (|2.11[ ). we can 
obtain the dimensions of the leading order perturbations 
for the DN fixed points. The basins of attraction cor- 
responding to the ND and DN fixed points thus show a 
mirror symmetry with respect to the line g c = g s in the 
interaction parameter space. 

Unlike the case of the junction of two quantum wires, 
where basins of attraction corresponding to the NN, DD, 
ND, and DN fixed points patch the whole space of the 
interaction parameters with some overlap, see Sec. IIV1 a 
section of the parameter space remains uncovered. This 
implies the existence of nontrivial fixed points. Notice 



that both ND and DN fixed points restore TRS even in 
the presence of the magnetic field. 



D. x+X+ and X-X- BC 

The most striking consequences of x+ anc f X- fixed 
point are that TRS is broken explicitly and the stabil- 
ity of the fixed points is determined both by the effective 
flux, cj), and the interaction parameters. For spinless elec- 
trons, the x+ fixed point can be stable when < (f> < tt 
while the X- can be- stable when — ir < 4> < 0- As a re- 
sult, the charge and spin current always flow together for 
spin 1/2 electrons, i.e. mixed chiral fixed points, x±X^> 
are never stable. Moreover, among all combinations of 
BCs corresponding to chiral BC, only x±X± fixed points 
are stable. One can conjecture that the criterion of sta- 
bility of the x± BC due to the flux for spin 1/2 electrons 
is the same as for spinless electrons. The boundary op- 
erators of leading order perturbations, such as backscat- 
tcring processes, have scaling dimension in terms of in- 
teraction parameters 



25c 



2g s 



5 2 



(2.12) 



and provide the only constraint for determining the sta- 
bility. The basin of attraction is thus common for both 
X±X± fixed points and is depicted in the Fig. [5] How- 
ever, which phase is preferred in the low energy limit is 
completely determined by the effective flux. 

The conductance tensors at the x±X± fixed points 
takes the Z% symmetric form in the Eq. 
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c(s) I 
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3 + s; 



[(3S jk -l)±g c(s) e jk \ (2.13) 



with the scalar conductance G 



c(s) 



antisymmetric component G 



cO) 



and the 



9c{s)G 
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The domi- 



nating hopping processes are schematically shown in the 
inset of the Fig. \5\ Notice that the conductance at 
the X±X± fixed points becomes perfect transmission of 
charge and spin from wire j to wire j ± 1 when g c ^ — > 1 . 



E. Asymmetric Boundary Conditions, DaDa BC 

We shall consider the simplest Z3 asymmetric BC, Da, 
corresponding to a situation where as shown in the inset 
of Fig [6l two of the wires are strongly coupled while the 
third one is decoupled from the rest of system. Note that 
the Da fixed point has been proven unstable in the Y- 
j unction system for spinless electrons 0- However, for 
spin-1/2 electrons, there is a window in which the D aD a 
fixed point is stable. Of course, it depends on the detailed 
structure of the hopping amplitudes to determine which 
wire will be decoupled. 
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FIG. 5: [Color Online] The painted area is the basin of at- 
traction of both x±X± fixed points. This area sits between 
and osculates the stable regions of the NN, DD, ND, and DN 
fixed points. The border shows the marginal line for scaling 
dimension A = 1 of Eq. p,12[ ). The pictorial representation 
of the tunneling processes associated with X+X+ fixed point 
is shown in the inset. An incoming electron from wire i will 
always divert to wire i + 1. The currents associated with the 
X_X- fixed point flow in the inverse direction. 
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FIG. 7: [Color Online] The proposed phase diagram - The 
colored lines indicate the marginal boundary of A = 1 for 
each fixed point upon the identification of the leading order 
perturbations. The painted areas shows the regions with only 
one stable fixed point while the unpainted regions represent 
those where two or three fixed points have an overlapping re- 
gion of stability. Notice that the DaDa fixed point is the sole 
stable one in four tiny gray areas surrounded by the unpainted 
regions. 
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FIG. 6: [Color Online] The painted areas are the basins of 
attraction for the DaDa fixed point. The overlap of the at- 
tractive basin for the DaDa fixed point and others is very 
complicated and will be illustrated in Fig. [7] A pictorial rep- 
resentation of the DaDa fixed point with the decoupled third 
wire is shown in the inset. Notice that there is a three fold 
degeneracy of this fixed point, i.e., the decoupled wire can be 
arbitrary, resulting from breaking the Z3 symmetry. 



As shown in Fig.0 the stable area of DaDa fixed point 
emerges from the borders of the attractive basins of other 
fixed points. Several boundary operators contribute to 
the leading order perturbations and need to be considered 
for determining the stability. Here, we will only present 
the basins of attraction corresponding to DaDa BC and 
postpone the discussion of scaling dimensions. 

The conductance at the asymmetric fixed points will 
not take the form in Eq. (|2.5p due to the broken Z3 sym- 
metry. The conductance tensor at the DaDa fixed point 
is given in the matrix representation 

e 2 ( 1 "I 0\ 

G c ^\ DaDa =2<7 C(S) - -1 1 . (2.14) 
11 \ / 

Despite an extra component representing the decoupled 
third wire, the conductance tensor is exactly the same 
as that at DD fixed point for junction of two quantum 
wires. 

We shall conclude this section by proposing a phase di- 
agram based on the results discussed above. The stability 
of phases at some area in the parameter space is uniquely 
defined because there exists only one stable fixed point. 
As a result, the painted areas in the Fig. [7| with red, 
green, blue, orange, purple and grey colors unambigu- 
ously correspond to the NN, DD, ND, DN, x±X±i 
and DaDa fixed points, respectively. When there are 
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more than one possible candidates for stable fixed points, 
determining the stability of phases becomes tricky and 
generically non-universal. Hence, the phases in the low 
energy limit are defined not only by the interaction pa- 
rameters, g c t s \ but also by the details of the device. For 
instance, the different hopping amplitudes of single elec- 
tron tunneling between wires may determine the final 
destination of the RG flow and the phase at the overlaps 
of DD and DaDa fixed points. However, determining 
the stable fixed points in these cases is beyond the scope 
of our methods and can only be conjectured. 

Let us consider the system along the SU S (2) invariant 
line, g s = 1: the NN fixed point is stable when g c < 
1, the chiral fixed points are stable when 1 < g c < 3, 
and the DN fixed point is stable when g c > 3. This 
result is expected due to the similarity with a system 
of spinless electrons where N fixed point is stable when 
g c < 1, the chiral fixed points are stable when 1 < g c < 3, 
and D fixed point is stable when g c > 3. Notice that 
four common tangential points of marginal lines when 
(9c,g s ) = (1, 1), (1,3), (3, 1), (3,3) remain marginal for 
all adjacent phases. 



III. MODEL AND BOSONIZATION 

In this section, we define the simplest model includ- 
ing the TLL effects and spin degree of freedom. As de- 
picted in Fig. [TJ we study three identical single channel 
quantum wires with spin 1/2 electrons joined with a ring 
enclosing a magnetic flux. We will ignore phonons and 
impurities, and assume that electron-electron interaction 
is short-ranged in the wires. The Euclidean action of a 
semi-infinite interacting single channel TTL with spin- 
1/2 electrons is described in terms of the independent 
charge and spin boson fields, 

3 

S = E / drdx{^[(d x(Pj , c f + ^{d T <p hC f] 



+ ^P^>) 2 + ^(^>) 2 ]}, (3.1) 

where the subscript j = 1,2,3 represent each wire and 
u c ( s ) is the sound velocity for the charge and spin degree 
of freedom respectively, or in terms of dual field 8 

p 1 

s = J2 / dTdx{-^-[{Mi,cf + -{d T e 3 , c f] 
+ ^[M, s f + ^{d T e j>s f]}. (3.2) 

Here, p c {s) an d # c (s) are phase fields which follow the 
canonical commutation relation, [<p c (s) [ x )i @c(s){ x ')] = 
— i®(x — x'), where 9(x) is the Heaviside step func- 
tion. Notice that g c = g 8 = 1 corresponds to the non- 
interacting point, and in the absence of a magnetic field 
and any spin-dependent interactions, we must take g s = 1 
in order to respect the underlying SU(2) symmetry [26j . 



The fields with spin up and down degrees of freedom 
can be represented as linear combinations of these charge 
and spin boson fields, 



tp-i + crtpsi 
(p at i = — — -j= , and 



<70„ 



V2 



(3.3) 



where the commutation relation between tp and 8 fields is 
still followed. It is convenient to introduce left and right 
mover representations 



(3.4) 



and to identify the electron fermion field in terms of the 
boson field 



,HR) 



(3.5) 



with Klein factors r]^ that satisfy the anti-commutation 
relation, {rj^,rj u } = 25^^ and commute with the boson 
fields. 

The effect of interactions at the boundary in the ac- 
tion Eq. (|3.1[) and Eq. (|3.2[) will be enter in the form of 
tunneling operators between wires. These boundary op- 
erators must conserve both charge and spin and therefore 
they are constrained to respect the corresponding U(l) 
symmetries. 



IV. DEBC 

The Delayed Evaluation of Boundary Condition 
(DEBC) method was introduced recently for determining 
the stability of the boundary conditions (fixed points) of 
the Y-junction for spinless electrons Q- In this section 
we will generalize this method to the case of spin-1/2 
electrons and use this method to determine the phase 
diagram of a junction of two quantum wires. 

Generically, fermionic operators can be represented in 
terms of the boson fields, p c (s) {x, t), and their conjugate 
fields, 6 c r s -\(x, t) up to a Klein factor in an infinite wire. 
For a semi-infinite wire, the relation between left and 
right moving fields leads to an analytical continuation 
4> R (~x,t) = cf> L (x,t). This is the familiar unfolded pic- 
ture where the right mover and the left mover are related 
by a particular choice of the boundary condition, namely 
the N BC, at x = 0. The DEBC method is an extension 
of the unfolded picture to different boundary conditions. 

Within the DEBC framework, an arbitrary boundary 
operator should be first represented in terms of the inde- 



pendent bulk boson fields 



MR) 



without specifying the 



boundary condition. Then the boundary conditions that 



must be imposed on 



MR) 



(t, x = 0) (cf>i and 8i) are deter- 



mined a posteriori. Because 4>^ R ^ are functions of x±t, 
imposing a particular boundary condition relates ip and 
8 fields in the bulk, and thus eliminates the redundancy 
of working with both tp and 9 fields in the semi-infinite 
wire. 
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The scaling dimension A of a boundary operator Ob, 
with two-point correlation 



(O, 



\ t -t'\ 



-2A 



(4.1) 



depends on the boundary conditions and can be used 
to determine whether a perturbation consisting of the 
boundary operators is relevant or not. When A > 1 the 
perturbation is irrelevant and when A = 1 the perturba- 
tion is marginal. Thus, a boundary condition is stable 
when all boundary operators either have scaling dimen- 
sion, A = 0, or are irrelevant A > 1. 

Now, let us introduce the generic representation of 
boundary operators and discuss how to obtain their scal- 
ing dimensions given a boundary condition. 



A. Boundary Operators and Scaling dimensions 

For applying the DEBC scheme, the left and right mov- 
ing fields appearing in the boundary operators initially 
are treated as independent regardless of the boundary 
condition. The boson representation of fermions is given 
in Eq. (j3~5)) 



*^(M = o) 



la-,! ° 



(4.2) 



All boundary operators are constructed by combining 
creation and annihilation operators of fermions in each 
wire. The Klein factors and extra phases produced by 
commuting the boson fields will not affect the scaling 
dimension of a boundary operator 



O, 



cxp 



(4.3) 



near stable fixed points. Here, n" ff uniquely defines the 
tunneling processes at the boundary, where i, er, and 
a represent the wire, the spin, and the chirality of the 
fermions, respectively. We shall refer hereafter to n° CT 
as the particle number vector. Since the boundary oper- 
ators are constructed from the full electron, the charge 



and spin degrees of freedom are coupled at the boundary. 
The total charge and spin conservation imply 



J2 < CT = 0, and ™l° = 



(4.4) 



respectively. For instance, a tunneling process where a 
up-spin right-mover at wire 1 scatters into a up-spin left- 
mover at wire 2 leads to n R ^ = 1 = — j with the sign 
convention followed from the bosonic representation of 
the fermions. 

It is convenient to introduce scaled bosonic fields 



<Pc(s) = ^9c(s)tPc(s) and 6> c ( s ) 



f C (s) 
\/9c(s) 



(4.5) 



such that the commutation relation between (p c [ s ) and 



is still followed. The action of the rescaled fields 



becomes independent of the interaction parameters, and 
their correlation functions are given by 



(0 c(s) (z,z)0 c(s) (O)) 



1 



lnlzl 



(4.6) 



In terms of the rescaled boson fields, the original left and 
right moving fields become 

V^i.tr = (cosha c (j) a i c + crcosha s </>° s ) 

+(sinha c (f>f c + asmha s <f>f s ), (4.7) 



where 



— , and 



l R - 

(7,1 



■ro 



(4.8) 



and cosha = (-^ + s/g)/2 and sinha = (-^ — ^fg)j2. 
Further, a is defined such that a = L when a = R and 
vice versa. Also, shorthanded notations, R — L and 
L = R, are used. 

In terms of the particle number vector and non- 
interacting bosonic fields, the bosonic argument of the 
boundary operators becomes 



J 



a / j 



{n a ia (cosha c 0° c + sinh a c ^ J +< (T a (cosha s $ s + sinha,; #,)}. 



(4.9) 



The non-trivial scaling behaviors of the boundary oper- are the fields 
ators are attributed to the mixing structure of the left 
and right movers in Eq. (|4.7p . 



^z,c(s) ' 



Eq. (|4~9)) becomes 



Let us introduce vectors 



TR(L) 



whose i components 



( v R ■ $ + v R ■ $ R ) + (v L c ■ & + • # ) 



(4.10) 
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where the vectors v R , \ ' are defined as 

c(s) 

(v R )i = /J(n^ CT cosha c + nf^ sinha c ) (4.11a) 

a 

("?)< = E CT ( n 5 cosha. + n^ sinha 8 )(4.11b) 

(T 

(«c)i = cosha c + < CT sinha c ) (4.11c) 

(T 

(«20< = E a ( n *> cosha.+n£, sinha s )(4.11d) 



BC corresponds to the total reflection fixed point and has 
7£;v = 1; The D BC corresponds to the perfect transport 
fixed point and has 



Tin = 



1 

1 



(4.15) 



Since the N and D BC can be imposed independently 
in the charge and spin sector, we can now discuss the 
contributions of the charge and spin degree of freedom 
to the scaling dimensions separately. 



Here 4> L = <\> R and <f> R = <fi L were used explicitly. Observe 
that the charge and spin degrees of freedom differ only 
by an extra a term. 

The boundary conditions can be generally identified as 



77?. _ -r> - 1 XL 



ft _ V + C|x=o, 



(4.12) 



where C is a constant vector and 1Z is a rotation matrix. 
The total charge and spin conservation imposes the N 
BC 



(4.13) 



on the center of mass mode, $o,c(«) = X)*Li 0i,c(» > 
which in turn constraints the rotation matrix 1Z. More- 
over, TZ has to be an orthogonal transformation to pre- 
serve the total length of the fields [2{|. The scaling di- 
mension A of the boundary operators for a boundary 
condition reads 



1 



A^« CT ) = -|ft 



(4.14) 



The physical processes responsible for pinning a 
boundary condition determine the rotation 1Z, for which 
Ak = for that specific particle number vector. The 
corresponding operators thus act as the identity opera- 
tor at this fixed point. Once the boundary condition is 
picked, the stability of the fixed point can be analyzed by 
evaluating the scaling dimensions of all other boundary 
operators. In particular, the fixed point will be stable if 
all these dimensions, for all other particle number vec- 
tors, satisfy A-]i{nf a ) > 1. 

To illustrate the method, here we apply the DEBC 
scheme to the case of a junction of two quantum wires 
for spin 1/2 electrons and find results in agreement with 
Ref. [i|[27L[2l. 



N BC in charge sector 



Focusing on the charge sector of Eq. (|4.14[) with 7Z C = 
1, the scaling dimension of the charge sector reads 



4 i 



(4.16) 



= — | (f^cr cosh a c + n1 a sinha c )e 



2 

e 2,c| , 



where we introduce ej iC as the basis vectors of v c and use 
cosha c = (-L + v /5c)/2 and sinha c = (^= - y^)/2 
explicitly in the second equality. Due to the charge con- 
servation, we can introduce 



5>k = -E«t 



(4.17) 



and the contribution to the scaling dimension reads 

1 



Wo 



-n 2 . 



(4.18) 



2. N BC in spin sector 

The contribution of the spin degree of freedom to the 
scaling dimension is very similar to the charge sector with 
a substitution of v s to v c in Eq. (|4.16p . Hence, the scaling 
dimension with N BC in the spin sector is given 



-\TZ s v R +v^\ 



(4.19) 

i2 



B. Boundary conditions of junction of two 
quantum wires for spin 1/2 electrons 

When requiring the spin and charge conservation and 
orthogonality, only two boundary conditions, Neumann 
and Dirichlet, are possible in the case of a junction of 
two quantum wires for spin 1/2 electrons. The 7Z matrix 
takes a particularly simple form in both cases. The N 



Due to the spin conservation, we define 



E 



E 



(4.20) 



with the relation n = s (mod2). The scaling dimension 
for the N BC on the spin degree of freedom reads 



1 

2g~ s 



(4.21) 
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3. D BC in charge sector the charge sector of the Eq. (|4.14p . we obtain 

The rotation matrix of D BC in Eq. (|4.15p implies the 
mixing of two quantum wires. Inserting the matrix into 



\U c v^ +v^\ = | E^ 71 ^ + n\ a ) cosha c + (n^ a + nf a ) sinhaj ei, c + + n^ a ) cosha c + {n\ a + nf^) sinha c ] e 2 , c 

(7 (7 

By using charge conservation, we introduce a new variable 

n = n + = J2 nl„ ©(+ ' «) = E«- + n V = ~ n ~ = ~ E n l° Q ^ 1 = ~ Efe + "a,«r)- ( 4 - 22 ) 

Here, when applied inside 9-function, the indices of wires I = 1, 2 become Z = +, — respectively, and the left and right 
mover indices become a = +, — respectively. In terms of the new variable, h, the scaling dimension is given by 

1 2 1 2 1 

Ad c = ^l^c^f + #c| = t|(VS^™) ei,c - Wfiln) e 2)C | = ( 4 ' 23 ) 

I 



4. D B.C in Spin Sector Eq. (|4.22j ). Since we have spin conservation rather than 

charge conservation, we can define a new variable 5 

Observe that the scaling dimensions attributed to the 
spin sector can be obtained by replacing v c with v s in 



= S+ = E °n?,o 0(+ I a) = E(™t + = - s " = - E e (- Z a ) = - E« + (4-24) 



I, a, a 



Also one can show the relation n = s (mod.2). The scaling dimension corresponding to the spin degree of freedom 
with the D BC is given by 



2 1 



r 



(4.25) 



C. Stability of the fixed points for junction of two 
quantum wires 



Since both the total charge and spin are conserved sep- 
arately, the boundary conditions can be imposed inde- 
pendently in the charge and spin sector. We follow the 
categorization in Sec. [II] by relating a fixed point to a 
B C B S BC, and explore its instability. The scaling di- 
mensions of boundary operators given a BC are obtained 
by adding the contributions from both the charge and 
spin degrees of freedom. In light of the discussion in this 
section, the possible scaling dimensions of boundary op- 



erators are given by 

And 
Adn ■ 



9c 9s 
g c n 2 + g s s 2 ] 

2 1 -2 

— n + g s s 

9 c 

g c n H s 

9s 



(4.26a) 
(4.26b) 
(4.26c) 

(4.26d) 



with the constraint n = s (mod 2) and h = s (mod 2). 
Notice that there is no constraint between n, s and be- 
tween h, s. 

Let us compute, as an example, the scaling dimension 
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of the following boundary operators given NN BC 

1. T B = ^t^i.t oc e^-^T+^i.r] (4.27a) 

2. T F = i, t a e iv/5[ -^T+^i, T l (4.27b) 

3. T = tyf^fj^^^ (4.27c) 

4. T = *^ T t *f i t *i, T *2.i- (4.27d) 

The case 1 corresponds to the backscattering process 
with = 1 and = —1, hence n = = s from 
Eq. p~T7)) and Eq. ljrfj)) . Notice that A% N = 0, and 
thus backscattering is the physical process that fixes this 
boundary condition. The case 2 corresponds to the for- 
ward scattering with = 1 and = — 1 such that 
n = 1 = s and Aj^ N = i(l/g c + l/g s )- The case 3 cor- 
responds to pair tunneling process with = 1 = nf\ 
and 7i^| = — 1 = r»2 j, One can show that n = 2 and 
s = and obtains A^ w = 2/g c . The last case corre- 
sponds to pair exchange tunneling process with n\ | = 1 , 

rii i = — 1) ^24. = 1 anc ^ n 2 t = — ^ such that n = and 
s = 2 and = 2/g s . The last three boundary op- 

erators describe the leading order perturbations for the 
NN fixed point. The basin of attraction of the NN fixed 
point is determined when all A2,3 ; 4 > 1. 

In principle, we need to know scaling dimensions of all 
boundary operators given a boundary condition in order 
to determine the basin of attraction for each fixed point. 
However, Eq. (|4.26p provides a compact way to determine 
the stability of the fixed points. Since the basin of attrac- 
tion is defined when the scaling dimensions of the pertur- 
bations have A > 1, we can obtain the possible scaling 
dimensions in an increasing order by inserting the small- 
est integers for the quantities n, s, n, and s with their 
respective parity constraints. Consequently, the stability 
of the fixed points can be completely determined from 
this minimum construction. The phase diagram deter- 
mined by DEBC for junction of two quantum wires for 
spin-1/2 electrons is depicted in Fig. [5] and agrees with 
the results in [H, [13, Hi . 

We conclude this section by reviewing the DEBC ap- 
proach. First, the boundary operators should be ex- 
pressed in terms of redundant fields ip and 9 or <j) R ( L > . 
Second, one has to identify the rotation matrix 7Z cor- 
responding to a given boundary condition. Finally, the 
scaling dimensions of boundary operators given a bound- 
ary condition are evaluated using the rotation matrix 7Z, 
and the basin of attraction of the fixed point can be de- 
termined by requiring all dimensions to be larger than 
one. 



V. JUNCTION OF THREE QUANTUM WIRES 
FOR SPIN-1/2 ELECTRONS-DEBC 

Following the DEBC scheme developed in the last sec- 
tion, we have to first identify the rotation matrix repre- 
sentations of boundary conditions for junctions of three 



4 









DD 

1 








1 


NN 








DN 



12 g c 



FIG. 8: Phase diagram of a junction of two quantum wires- 
The shaded area is common to the NN and DD stable regions. 
There should be an unstable intermediary fixed point in this 
range of g c{s y 

quantum wires for spin 1/2 electrons. The conservation 
of the total charge and spin at the boundary can be im- 
plemented by imposing the N BC on the zero mode (total 
charge or spin mode) as shown in Eq. (|4.13[) . Also, the 7Z 
has to be orthogonal to preserve the length of the field 
vector. Observe that any rotation around the unit vec- 
tor, uq = ^=(1, 1, 1), will leave the total charge and spin 
mode invariant, hence naively 7Z can be an arbitrary rota- 
tion around uq- However, there are some boundary con- 
ditions corresponding to discrete transformations, where 
1Z cannot be totally classified by a rotation. Moreover, 
as discussed in the Appendix [Al the conductance of fixed 
points can be directly related to the rotation matrices, 
Eq.flllUs 

e 2 

Gij,c(s) = 29c(s)-r($ij - TZji)- (5-1) 

We will compute the conductance by the construction of 
the rotation matrices associated with different boundary 
conditions. 



A. Boundary Conditions 

We will construct in this subsection the transforma- 
tions 7Z corresponding to different boundary conditions. 
Generically, any rotation matrix satisfying the constraint 
of conservation and orthogonality can be a candidate for 
a physical boundary condition. However, only some of 
them correspond to stable fixed points. Below, we will 
calculate the TZ matrices for physically motivated bound- 
ary conditions. 
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1. Neumann boundary condition 
The Neumann boundary condition, 

J^(z)U=o = <=> 9\ x=0 = , (5.2) 

corresponds to a fixed point where the three wires are 
decoupled from each other. The relation 6\ x= o = 
can be represented in terms of the left and right mov- 
ing fields, <j) R = (f> L \ x =o- The backscattering processes 
thus yield zero scaling dimension and dominate the low 
energy physics. The rotation matrix of the N BC is given 
K N = 1. 

By using Eq. (|5.ip and the rotation matrix of N BC, 
the conductance is 

Gt ) N=^c( s )^(S ij -S ji )=0 (5.3) 



3. Chiral boundary condition 

It is trickier to obtain the transformation correspond- 
ing to the chiral fixed points since the rotation matrix 
depends on the coupling constants. However, the pinning 
processes with A = for the chiral fixed points provide 
the conditions for constructing the rotation matrix. 

Let us consider the charge sector of the chiral-like tun- 
neling processes: A left mover at wire i scatters to a right 
mover at wire i+p, i.e. nf a — 1 and nf~ +p a = — 1 for all 
i. The scaling dimension of these boundary operators for 
arbitrary rotation matrix can be written using Eq. (|4.14l) 

A = \1Z (sinha c — cosha c ii +p ) 

+ (cosha c ii — sinha c £i+ p ) | 2 . (5-8) 

By using the fact that A x+ = for all chiral boundary 
operators with p = 1 and A x _ = for those with p = 2 
and after some algebra, one obtains the rotation matrices 



2. Dirichlect B. C-Andreev reflection fixed point 

Since the D BC, <p(x)\ x =o = C , cannot be imposed on 
all boson fields due to the N BC imposed on the center of 
mass mode, it can be only imposed on two independent 
fields 

$\x) = -L(^-^ 2 )(x)U =0 = 0, (5.4a) 
$ 2 (x) = -±=( (pi + 92 -2 i p 3 )(x)\ x = = 0, (5.4b) 

where $ 1 and <E> 2 are orthogonal to the center of mass 
mode. Hence, the D BC is given by the condition 



(5.5) 



This transformation can be described by a rotation 
around So by n, and is given by 



1 


2 


2 


2 3 


3 1 






2 3 


\ 


3 


3 


3 



Then the conductance of the D BC is 

,2 1 



(5.6a) 



(5.6b) 



G?A = 2 5c (a)^(*ij - 3(2 - 3^-)) (5.7a) 



2^^(3^-l), 



(5.7b) 



with G 



c(s) _ 8g c(s) 



and 





2(l+9c) 


2(l-3c) 


3+9? 

2(l-5e) 


3+9? 
-1+9 2 . 


3+9? 
2(l+9e) 


3+9? 

2(l+5c) 


3+9? 

2(l-ffo) 


3+g? 
-1+9 2 


3+9? 


3+g? 


3+9? 


-l + Sc 


2(l-flc) 


2(l+9c) 


3+9? 

2(l+ffc) 


3+9? 


3+9? 
2(l-9c) 


3+9? 

2(l-9c) 


3+g? 

2(X+ flc ) 


3+9? 
-l+9 r 2 


3+9? 


3+9? 


3+9? 



(5.9a) 



(5.9b) 



or in the tensor form 
1 



3 + <r 



[(-3 + 5 ^+2(l± 5 ey)], (5.10) 



for x± BC respectively. Similarly, the 7?. matrix for the 
spin sector will be in the same form as Eq. (|5.9p with 
a simple substitution, g c — > g s . For the case of more 
than three wires, there may exist more than two "chi- 
ral" boundary conditions. In this paper, we will restrict 
oursclf in the case of three quantum wires. 
Using Eq. (|5.1[) . the conductance is given by 



G 



:(s) 



h 



3)Si,j + 2{l±g c{s) e jti ) 



h 3 + g 



■[(3S ld - 1) ±g c ( s y 



(5.11) 



4- Asymmetric fixed points, Da 

There are boundary conditions that cannot be repre- 
sented by a simple rotation around the Uq; instead, they 
correspond to improper rotations. The only important 
fixed point that falls into this category corresponds to the 
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asymmetric boundary condition (Da BC), in which one 
of the wires is totally decoupled from the system. The 
matrix representation of this boundary condition can be 
constructed by imposing the N BC on the third wire and 
the D BC on the first and second wires. ( Here, one can 
choose to decouple any of the three wires. Hence, there 
are three distinct D\ BC, where i = 1,2,3 represents the 
decoupled wire. However, they are identical up to a Z3 
transformation 
given by 



The matrix representation of D\ BC is 



cycle: 

= VfJ^tloce-i^e^*, (5.15a) 



2.<r 



2,<x|0 



,R U 



:>>.a 



2, a- ^ 3,cr I q 



OC e 



(5.15b) 
-, (5.15c) 



and also classified by the non-zero elements of particle 
number vector for S~ a : n\ a — 1 and nf_ l a = — 1. 
• Backscattcring: 




(5.12) 



According to the Eq. (|5.1[) . the conductance tensor of the 



D\ BC becomes 



G 



c{s) 



2ffc(s)- 








(5.13) 



Now, we can impose these boundary conditions sep- 
arately on charge and spin sectors and determine the 
stability of each fixed point by computing the scaling 
dimensions of the boundary operators. 



B. The boundary operators 



oB 

°1,<T 


= *v<« 





OC 


e~ 






oB 
°2,a 


= 





OC 








oB 
D 3,a 







OC 


e~ 







(5.16a) 
(5.16b) 
(5.16c) 



Again, we should identify the representation of the par- 
ticle number vector for Sf a : n\ a = 1 and nf- a = — 1. 
• LL — RR combinations: 



oL 



J 3,a 



nit 
1,0 

oR 
2,0 

oR 
d3,<r 



l,cr 3,0 
,R t, T ,i?, 



(5.17a) 
(5.17b) 
(5.17c) 



** e iV2<l>? ia ( 5 .17d) 



K0 K0I « e- iV ^e^ % (5.i 7 f) 



Now, we have to write the boundary operators in terms 
of redundant bosonic degrees of freedom. In principle, we 
need to consider all tunneling operators between three 
wires at x = 0. However, since more particle processes 
tend to be less relevant, we will merely consider single and 
two particle tunneling processes. The scaling dimensions 
of these operators given the BC will be evaluated later 
and provide measures of relevant or irrelevant perturba- 
tions. A more systematic method to compute all scaling 
dimensions of boundary operators, which relics on con- 
formal symmetry of the system, will be introduced in the 
next section. From the full spectrum of the dimensions, 
one will conclude that the leading order perturbations 
indeed come from the few particle tunneling processes. 

The single electron tunneling processes can be classi- 
fied into four groups and their Hermitian conjugates: 

• + cycle: 



Then, we can identify the representation of the particle 
number vector for Sf^ : = 1 and = — 1. 

Here, S indicates "single" particle processes. Notice 
that there can be no spin flips in single particle processes. 
However, the scaling dimensions of multi-particle tun- 
neling processes depend on the spin degree of freedom. 
The multi-particle operators can be constructed from the 
combinations of single particle boundary operators and 
their hermitian conjugates. Here, we list some of the two 
particle processes that will be useful for identifying the 
leading order perturbations. 

• Pair Tunneling in + cycle: 



(5.18a) 
(5.18b) 
(5.18c) 



PT+ 




PT+ 


— o 2 ,T°2,i 


PT+ 


- 9+ S"+ 





= *VK 




= ^3 J ^2. 




= *iV*3. 



(5.14a) 



K g-tVS^gtVS^ ( 5 14b ) 



(5.14c) 



These boundary operators can be categorized by the non- 
zero elements of particle number vectors n° CT for Sf a : 
-1. 



n Z,a = 1 and n ?+l,a 



Then, we can identify the representation of the particle 



number vector for PT„ 

1.. 7 - 1 - 
• Pair Tunnc 



+ . n L 



= nf, = 1 and 



in — 


cycle: 


PTf 




PT-2 


= s 2 ^s 2A 


PT 3 - 


= 83,^3,1 



R 

(5.19a) 
(5.19b) 
(5.19c) 
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In terms of the particle number vector, we obtain for 



pt; 



1 and nf_ 1 



R _ 



-1. 



• Pair Tunneling in LL-RR combinations with net spin: 



ryrp qLR, qL qR 

r 1 °\,o — °l,a°l,<7 

ryrp qLR qL qR 

r 1 °2,er — D 2,a D 2,a 

td^tqLR. qL qR 

rl °3,<t — a 3,o a 3,a 



(5.20a) 
(5.20b) 
(5.20c) 



In terms of the particle number vector, we obtain for 



PTS, 



LR. L 



n 



R _ 



+1,<T 



= -1. 



• Pair Tunneling in LL-RR combinations without net 
spin: 



TjrpLR, qL qR 

rJ L(r — l,a°l,-a: 

TjrpLR qL qR 

rJ 2,tr — D 2,a>~>2,-(ji 

jyrpLR qL qR 

ri 3,tr — '- , 3,cr' : '3,-o-- 



(5.21a) 
(5.21b) 
(5.21c) 



One can read the the particle number vector off for 
PTt*: n^ = nl_ a = 1 and nf +1>a = nf +h _ a = -1. 
• Pair Backscattering in the same wire: 



PB\ = Si^S^ 
PB 2 = 

PB 3 = S 31 S 3l 



(5.22a) 
(5.22b) 
(5.22c) 



The representations of non-zero elements of particle num- 



1 and nf| 



ber vector become for PBi 

• Pair Backscattering in the different wires without net 
spin: 



PB 12 , a 


oB qB 
— l,£r 2,-CT 


(5.23a) 


PB 23 . a 


qB qB 

— °2,<r°3,-CT 


(5.23b) 


PB 31 . a 


qB qB 
— °3,<J°l,-a 


(5.23c) 



The non-zero elements of the particle number vector are 
given for PB io , a : nf a = nf _ a = 1 and nf a = nf_ a = 
-1. 

• Pair Backscattering in the different wires with net 
spin: 

PBS l2 , a = Sl a Sl a (5.24a) 
PBS 23 ,a = Sl„Sl a (5.24b) 
PBS 3ha = S*S* (5.24c) 



The non-zero elements of the particle number vector are 
given for PBS, 



Pair Exchange processes: 



n la = 1 and n f,a 



n R = -1 



PE h<7 


= ^l,a^2-a 


(5.25a) 


PE 2 , a 


= ^2,a^3,-a 


(5.25b) 


PE 3 , a 


= ^3,a^l,-a 


(5.25c) 



The non-zero elements of the particle number vector are 



given for PEi^: n. 
nf „ = -1. 



L 

La 



1 and nf +1 



Particle Hole Pair Tunneling in + cycle: 

(5.26a) 
(5.26b) 
(5.26c) 

The non-zero elements of the particle number vector are 
given for PH^: nf^ = j_ = 1 and nfy 
-1. 

• Particle Hole Pair Tunneling in — cycle: 



PH+ 




PH+ 


— D 2,]°2,i 


PH+ 


— °3,T 3,1 



PH 1 = S l1 -S l ^ 
PH 2 = S 2 ^S 2 ^ 
PH 3 = S 31 -S 3i 



(5.27a) 
(5.27b) 
(5.27c) 



The non-zero elements of the particle number vector are 
given for PH~: nf T = rtfLi x = 1 and nf l = nf_ x T = 
-1. 

• Particle Hole Exchange processes: 

(5.28a) 
(5.28b) 
(5.28c) 

Here, the non-zero elements of the particle number vector 



PHE ha 


= ^l,a^2, 


PHE 2 , a 


= ^2,a^3, 


PHE 3 , a 


= ^3,a^l, 



are given for PHEj „ : nf „ 



1 and nf"„ 



H+l,-a 



= -1. 



Using the particle number vectors of each boundary 
operator, the dimensions of each operator given the BC 
can be evaluated using Eq. (|4.14p . Hence, the stability of 
each BC can be determined by considering the relevant 
or irrelevant perturbations induced by the boundary op- 
erators. 



C. Stability of fixed points 

This section is organized in terms of the different 
boundary conditions. The boundary conditions are im- 
posed separately on the charge and spin degrees of free- 
dom. We will evaluate the scaling dimensions of the 
boundary operators corresponding to the leading order 
perturbations. 



1. NN BC 

The NN BC corresponds to the totally reflective fixed 
point in both charge and spin degrees of freedom. Hence, 
all backscattering processes, Eqs. (j5~T! IST221 15T231 [5T24]) 
and their hcrmitian conjugates will have zero scaling di- 
mensions. Now we identify the boundary operators asso- 
ciated with the leading order perturbations. 
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All other single particle processes contribute to the 
leading order perturbations with scaling dimension 



1 



1 



A s - i 

2g c 2g s 



(5.29) 



The pair tunneling operators in the + and — cycle, 
Eqs. (|5.18p and (|5.19p . the pair tunneling in LL-RR oper- 
ators without net spin degree of freedom in Eq. (|5.21[) and 
the particle- hole exchange operators in Eq. (|5 . 28[) are also 
leading order perturbations with the scaling dimension, 



. PT* 



, PT 



LR ■ PHE _ ^ 



*NN 



= A 



NN 



9c 



(5.30) 



Finally, the pair exchange operators in Eq. (|5.25[) have 
scaling dimension 



a pe 

^NN 



(5.31) 



and are therefore leading order perturbations. 

Now, we can identify the stable region when all leading 
order perturbations are irrelevant, A > 1, and obtain the 
basin of attraction shown in Fig. [2] 



2. DD BC 

The DD BC corresponds to the Andreev reflection 
fixed point in which all pair tunneling LL-RR operators 
with net spin, Eq. (|5.20p . have zero scaling dimension. 
Note that the dominant processes are multi-particle op- 
erators. 

First, the single particle processes in the + and — cy- 
cles in Eqs. (|5.14p and (|5.15[) are leading order operators 
and have the scaling dimension 



DD 



1 



(9c + 9s)- 



(5.32) 



In addition, the pair tunneling operators in ± cycles 
in Eos. (|5. 181 I5.19[) . and the pair exchange operators in 
Eq. (|5.25[) are also leading order perturbations with scal- 
ing dimension 



A 



PT ± 
DD 



:9c 



(5.33) 



Moreover, the particle-hole pair tunneling operators in 
Eqs. (|5.261 I5.2T[) and the particle hole exchange operators 
in Eq. (|5.28p have scaling dimension 



A PHE, _ A PH? 
^DD — ^DD 



= 3 9s 



(5.34) 



Thus, the basin of the attraction corresponding to the 
DD fixed point can be obtained by requiring A > 1 and 
is shown in Fig. [3J Unlike the case in junction of two 
quantum wires, the stable region of NN BC and DD BC 
do not overlap. 



3. ND BC 

The ND BC corresponds to a fixed point which de- 
scribes a charge insulator but spin conductor. The pair 
exchange operators in Eq. (|5.25p fix the boundary condi- 

PE 

tion with A N £'" = 0. Intuitively, the pair exchange pro- 
cesses induce a pure spin current depicted in the inset 
of Fig. 0] Furthermore, the processes of pair backscat- 
tcring in the same wire Eq. (|5.22[) arc pinned as well and 
indicate that there is no charge current. 

The leading order perturbations are attributed to the 
operators for the single particle tunneling in the ± cycles, 
Eqs. (|5.141 15.15p . with the scaling dimension 



1 

25c 



A s± - — -I- SI 
^ND — ~ r — 



9s 
6 



(5.35) 



The single particle backscattcring processes in Eq. (J5.10D , 
the pair backscattering processes in Eq. (|5.24p and the 
particle hole pair tunneling in Eqs. (|5.26[ I5.27P are also 
leading order perturbations with 



A S B _ aPBS 



3 



(5.36) 



Moreover, the pair tunneling operators in Eqs. (|5. 18115. 191 
I5.20p are also leading order perturbations with scaling 
dimension 



PT 



_ \PTS L 
— ^ND 



2 

9c 



(5.37) 



Again, we can identify the basin of attraction for the 
ND fixed point by requiring all A > 1 as shown in Fig. [4] 



4. DN BC 

The DN BC is the counterpart of the ND BC and cor- 
responds to a fixed point where the system becomes a 
charge conductor but spin insulator. Here, the particle 
hole exchange operators in Eq. (|5.28p fix the boundary 
condition with A^^ 1 '" = 0. As illustrated in the inset 
of Fig. [4] there is no spin current in these tunneling pro- 
cesses. Moreover, the pair tunneling operators in LL-RR 
channel without net spin Eq. (|5.2ip are also pinned at DN 
fixed point since they also represent the processes with 
only net charge current. 

The leading order perturbations are attributed to the 
single particle processes in the ± cycles in Eqs. (|5.141 
I5.15P with scaling dimension 



A s± 

^DN 



9c 
6 



1 



(5.38) 



The backscattcring in Eq. (|5.16p . the pair tunneling op- 
erators both in ± cycles in Eqs. (|5.181 I5.19P and the 
pair backscattering processes involving different wires in 
Eqs. (|5.231 15.24|) are leading order perturbations with 



^DN 



^DN 



P Bij 
DN 



PBSi 
^DN 



2gc 

3 ' 



(5.39) 
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Finally, the pair tunneling operators in LL-RR channel 
with the net spin in Eq. (|5.20[) and the particle-hole pair 
tunneling in both ± cycles in Eqs. (|5.26l 1 5 . 2 T[) represent 
the other set of leading order perturbations with scaling 
dimension 



A PTS 



PH ± 



2 

9s 



(5.40) 



Again, the basin of attraction can be found for the DN 
BC by requiring all A > 1 and is shown in Fig. [4j 



5- X+X+ and X-X- BC 

The X+X+ BC describes a fixed point in which both 
charge and spin current have a preferred flow direction 
1 — > 2, 2 — > 3 and 3 — > 1. Thus, at the X+X+ fixed 
point, the processes involving current flows in these par- 
ticular directions will fix the boundary condition. For 
instance, the pure charge current processes in the + cy- 
cle, Eq. (|5.14[) and Eq. (|5.18[) . and the pure spin current 
processes in the + cycle, Eq. ([5.26[) . have zero scaling di- 
mension. 

Notice that the single particle tunneling operators in 
the + cycle have zero scaling dimension regardless of 
the spin degree of freedom. However, boundary oper- 
ators representing pure spin current, constructed from 
the particle hole pair tunneling, have zero scaling dimen- 
sion only when the spin current is in the + cycle, i.e. 
A PH 7^ 0. Similarly, operators corresponding to the 
pure charge current, constructed from pair tunneling pro- 
cesses, have zero scaling dimension only when the charge 
current flows in the + cycle, i.e. A PT ^ 0. 

One can confirm that all single particle processes ex- 
cept the + cycle and some multi-particle processes pro- 
vide leading order perturbations with scaling dimension 



A 



2<? c 



2<? s 



x+x 4 



3 + gi 



(5.41) 



The scenario of X-X- BC is very similar to the case 
of X+X+ with relative changes from the + cycle to the 
— cycle. The leading order perturbations have the same 
scaling dimension as that of the X+X+ BC. Thus, both 
fixed points have exactly the same basin of attraction 
in the coupling constant space. The stability of the two 
fixed points is determined by the direction of the mag- 
netic flux threaded through the ring. We then plot the 
stable region of the fixed points in Fig. [5] 



6. D A D A B.C 

Without lose of generality, we will only consider in this 
section the boundary condition where the third wire is de- 
coupled, D A , both in charge and spin degrees of freedom. 
The situations with the first or second wire decoupled are 
similar. 



In the presence of the Z% symmetry, the operators in 
the same group listed in the previous section have the 
same scaling dimension. However, the DaDa fixed point 
breaks the Z% symmetry, hence the operators in the same 
group may have different dimensions. For instance, the 
operators in Eq. (|5.14al I5.15ap , have zero scaling dimen- 
sion and fix the BC, while the rest of the operators in the 
± cycle are leading order perturbations with the scaling 
dimension 



A = 



9c 



3 + ff s 2 
8<? s 



(5.42) 



In addition, the single particle backscattering in the third 
wire in Eq. (|5.16c[) has zero scaling dimension at the 
D A D\ fixed point while the backscattering at the first or 
second wire becomes another leading order perturbation 



A = i( 5c 



9s 



(5.43) 



Some of the operators associated with the leading order 
perturbations are listed below: 



Operators 


Dimensions 


pr cLR 
rl D 2(3),o- 


HSc 1 + 97 1 ) 


PT 2(3)> PT 2i3),a' PHE 2(3),a 


3 , g 3 

2 Sc 2 


PTf*, PB 12 , a 


2.9s 


PB m 


2ff c 




9a 1 3 

2 X 2g s 



(5.44) 



Thus, we identify the basin of attraction for DaDa BC 
and plot it in Fig. [S] 



7. NDa and DaN BC- A demonstration of the unstable 
fixed points 

In principle one can arbitrarily combine different BC in 
charge and spin sector to obtain the new boundary con- 
ditions. However, the rest of them are unstable against 
perturbations. We discuss here two unstable fixed points, 
NDa and DaN BC, and show explicitly that there are 
always leading order perturbations with scaling dimen- 
sions smaller than one in any region of the interaction 
parameter space. Considering the case where the third 
wire is decoupled from the ring, the backscattering in 
the third wire fix the BC in both ND\ and D\N fixed 
points. The operators PH^ and PH^ , corresponding to 
the processes with pure spin current between the first and 
second wire, have zero scaling dimensions at the ND\ 
BC while the operators PT^~ and PT^ , corresponding to 
the processes with pure charge current, have zero dimen- 
sion at the D A N BC. Hence, we conclude that ND\ BC 
corresponds to a fixed point with pure spin current be- 
tween first and second wire, while D A N BC corresponds 
to a fixed point with pure charge current between them. 
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We list below some operators which are crucial for de- 
termining the stability of the fixed points with their scal- 
ing dimensions, 



(5.45) 
(5.46) 



oB 

S 1(2),<t 


9s 




3 


ND\ 


2 ' 


~ 2^1 


oB 

D\N 


9c 
2 ' 


r>rpL R 

A- r - i 2(3),<7 


3 

" 2fc' 



Observe that A JS'' 

NDi 



and A ND , 



cannot be larger than 

s B . 

one at the same point in the parameter space; A^ 

and A D3 2 ^''" cannot either. Hence, we can concludes 

that NL>a and DaN fixed points are unstable against 
the perturbations. 



the standard free boson. The twisted structure affects the 
possible boundary states, the scaling dimensions of the 
boundary operators and the stability of the fixed points. 
In this section, we will first review the mode expansion 



of bosonic fields <p c r a \ and 



Ms) 



and derive the twisted 



structure of a TLL quantum wire for spin-1/2 electrons. 
A simplified version of the BCFT, exclusively for non- 
resonante tunneling Q, will be implemented by project- 
ing out the degree of freedom corresponding to the to- 
tal charge and spin. As a demonstration, we compute 
the partition function and obtain the scaling dimensions 
of the boundary operators corresponding to a boundary 
state for the case of a junction of two quantum wires. 



Review of mode expansions of boson fields and 
derivation of twisted structure 



D. Summary of DEBC 

We have demonstrated in this section how to imple- 
ment the DEBC method for obtaining the scaling dimen- 
sions of the boundary operators and determining the sta- 
bility of different fixed points. We find that the DEBC 
method provides a simple way to examine the junction 
systems and to determine the phase diagram. However, 
the main drawback of this method is that one has to 
determine the scaling dimensions of "all" boundary op- 
erators in each given BC. As mentioned in this section, 
we have to rely on the conjecture that operators involving 
more particles are less relevant to simplify the computa- 
tion. Hence, in the next two section, we shall provide 
the confirmation of our results through the approach of 
boundary conformal field theory where the full spectrum 
of the scaling dimensions of the boundary operators can 
be identified explicitly. 

VI. REVIEW OF BCFT 

The application of boundary conformal field theory 
(BCFT) to the analysis of critical phenomena was devel- 
oped by Cardy and widely applied to the Kondo problem 
and one dimensional problems [10, [U, [HJ . In the BCFT, 
comformally invariant boundary conditions are formu- 
lated in terms of boundary states. Here, the boundary 
states belong to the Hilbert space of the theory with the 
periodic boundary condition in the space direction. Con- 
sidering the Fermi statistics of the electrons, the struc- 
ture of the Hilbert space is twisted compared to that of 



The BCFT of free bosons can be applied to the system 
of interacting electrons via bosonization. However, due 
to the Fermi statistics of the electrons, there are various 
differences compared to the standard compactified boson 
field theory. In this section, the effect of the Fermi statis- 
tics will be implemented in terms of a twisted structure 
in the Hilbert space of the free boson theory. 

Since the Matsubara formalism for fermions implies 
anti-periodic boundary condition in the imaginary time, 
the space direction will also have anti-periodic boundary 
condition after the modular transformation. The bound- 
ary condition should be imposed independently on left 
and right movers, 

4>R(L),*{t, X) = -l/jR(L) !<T (t, X + /3) . (6.1) 

In the case of free fermions, g = 1, the boundary condi- 
tion becomes 

exp[i^ fl( £ )j0 .(t, 0)] = - exp[iV2^ R{L)tCT (t, 0)] . (6.2) 

Since [4>r(l)i 4>r(L)} 7^ 0, it is not clear how to determine 
the boundary condition on the boson field from Eq. (|6.2[) . 
We refer to Ref. 0] for the proof that the anti-periodic 
BC on fermions leads to the periodic BC on bosons. Pre- 
cisely, the right and left moving boson field, (j)R(L) are 
periodic variables 

4>R(L),v{t, % + 0) ~ 4>R(L),a(t, x) + V2irn, (6.3) 

where (n € Z). The mode expansion of the boson fields 
compactified on a circle of radius l/v2 becomes 
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Mx-) = <$„ + ^Q?x. + ±= £ -U<<. e- i2 F*- + h.c] (6.4a) 

where the momentum variable conjugating to <P(y^ , have eigenvalues \/2 times an integer, cin^ are 
oscillator modes, and x± = t ± x. 

Following the standard bosonization scheme, the independent charge and spin modes are defined as 

<t>R{L),c = ^(0R(i),T +<t>R{L),i) (6.5a) 

4>R{L),s = -^{4>R{L)A - <t>R{L),[) ■ ( 6 ' 5b ) 

Hence, the mode expansions of the boson fields corresponding to the charge and spin degrees of freedom become, 



oo 



6k,c( S ) = 0O lC(s) + ^Qc(s)t- + ^= E ^KcW e l ^ X ~ + h - C -1 ( 6 ' 6a ) 



71=1 

with the definitions 



(*) = <P0,c(s) + -^Qc(s)t+ + E ^ a n,c(s) e * 2 '" X+ + h A > ( 6 - 6b ) 



£ R(L) _ J_ ( .R{L) .R(L), R{L) _ J_, R{L) , R{L), , AR(i) _ AR(L) , Afl(i) , fi fi x 

^0,c(a) _ T/f^O.T ± ^0,i J' Vc(s) ~ TTf 1 -",! ± °X1 )' &na *°t c(s) ~ ^| ± ^J. • 



Notice that the eigenvalues of the operators Q^^ follow Now, let us define the new boson fields in the new basis, 
a special relation 

fc{s) = 4>R.c(s) + 4>L.c(s), 9c(s) = 4>L.c(s) ~ 4>R.c(s)- (6-8) 

Q R{L) =Q R(L) ( mod2 ), (6.7) 

Following Eq. (|6.6a|) and Eq. (|6.6b[) . the mode expansion 
named gluing conditions. of the new felds is given by 

I 



1 1 

Vc(s) = 0o,c(.) + + Qc{s)x] + 7| E ^Kcwe - ^*" + «t( s )e"'T I+ + h.c] (6.9a) 

* M — 1 * 

7T ' 1 °° 1 2 2 

O (.) = 0o,c(s) + -g[Qc(.)* + Qc(s)^] + 7| E ^[OM.)^^ - <c(,)e" !;2fia: - + h.c] , (6.9b) 



where 

0O,c(s) = 0O,c(s) + ^c(a)) ^0,c(«) = 0O,c(s) _ ^c(s) ) Qe(s) = Qc( s ) + <5f( s ) , <2c(s) = <2c( s ) ~ ■ (6.9c) 



A new gluing condition can be derived from Eq. ()6.7j) Hence, the periodic boundary condition of the fields <^ c ( s ) 



and 9 c r s \ are of the for 



111 



Qc = Q s = Qc = Qs (mod2), (6.10a) 

Qc + Qs + Qc + Qs = (mod4). (6.10b) </j c(s ) - p c ( a ) + 7rQ c(s) , c(s) - c(s) + 7r(5 c(s) , (6.11) 
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accompanied by the gluing condition in Eq. (|6.10[) . We 
refer to this as a "twisted structure", which reflects the 
fact that the bosons arise from the bosonized fermions. 
We have only discussed the mode expansion of the free 



fermions, g = 1, so far. However, similar results can be 
carried over to the case of interacting fermions with arbi- 
trary g. The mode expansions of yw s ) and 6 c u\ become 



7T r 1 



c(s) 



lc(s)t + Qc(s)X\ + 



E ^{< c{s) e- l ^ x - + < cW e- i2 ? I+ + h.c] (6.12a) 

n— 1 V 



7 C {» 



Vc( S ) + ^ [<5 C («)< + 5c( S )Q C ( S )^] + y ^ E 7^ 

V ^ n=l * 



a^e-^^-a^e-^^+^c] , (6.12b) 



following the same periodic boundary conditions in 
Eq. (|6.11[) and gluing conditions in Eq. (|6.10[) . 



B. Junction of two quantum wires for spin-1/2 
electrons 



In this subsection, we analyze the stability of the junc- 
tion of two quantum wires for spin-1/2 electrons by com- 
puting the partition function and the scaling dimensions 



of the boundary operators provided that charge and spin 

are conserved and there is no resonant tunneling. This 

method has been pursued by Wong and Affleck [28| for 
a junction of two quantum wires; here, we follow closely 
their approach, and only differ in that we make use of the 
conservation laws at the very beginning, thus changing 
the gluing conditions and simplifying the computation of 
scaling dimensions obtained from the expansion of the 
partition function. 



J 



1. finite-size spectrum, and bulk operators 

Let us consider two sets of independent boson fields Ve(«) an< ^ ^c(s)' wnere i = 1, 2 labels the wires 1 and 2. Since 
total charge and spin are conserved, it is convenient to work in an alternative basis, 



4> 



o 



1 



(6.13) 



where the total charge and spin modes are explicit. Using the Eq. (|6.12a[) and Eq. (|6.12b[) . the mode expansions of 
and < &c(s) become 



= M , s 

c(s) °' c(s) ' V2/3 l g. 



[^ + 0°( S )*] + 



y -L[ a ^° ( ,e^ x - + a Lfi , ,e- l2 -^ x + + h.c] (6.14a) 

,f^ v n,c(s) n,c(s) J V ) 



®c(s) = $0,c(») 



7T r Q 



c(s) 



0) 



t + Qc(s)X] + 



1 y J_ r ^ 

y/ 2 9c(s) £j 



a n,c{s) 



e b 



a n,c(s) e 



l-g-X + 



h.c] , (6.14b) 



where 



^R/l,o _ _~_/n/u,±, ji/u,£\ a njl " = - (a n ' llJ ^ af 1 '^'*) 

n,c(s) /o 7l,c(s) ro,c(s) •'' n,c(s) /o H,c(s) n,c(s)' : 



R/L,l R/L,2 



R/L 



1 / R/L,2s 



V2 



V2 K 



(6.14c) 
(6.14d) 



and the winding operators follow 

Qc{s) = Qc(s) + Qc(s)' Qc(s) = Q c {s) ~ Qe{ S y Qc(s) = Ql{s) + Ql(s)> Qc(s) = QI( s ) - QI( s )- (6.14e) 
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The eigenvalues of these winding operators are integer and follow gluing conditions, derived from Eq. (l6.10p : 



1. Q° c = Q° s = Q° c = = Q c = 

2. Q° + Q° c + Q° s + Q° s = 

3. Q° + Q° C +Q S + Q C = 

4. Ql + Q° c + Q s + Qc = 
5- Q° c + Q° s + Q° c + Q° s + Q 



Q s + Qc + Qs = 



(mod2) 
(mod4) 
(mod4) 
(mod4) 
(mod8) 



(6.15) 



Using the Lagrangian density 



3 = 1,2 



-in 



the Hamiltonian of the two quantum wires can be written in terms of the winding and number operators as follows 



2n 
J 



(Qg) 2 , g c (Qg) 2 , (Q c ) 2 , g c (Q c f , (Qg) 2 , g s (Q° s ) 2 , {Q s f , g s (Q s ) 
16 g c 

OO 



16 



•il,0 



163c 

-i,0 



16 16g s 16 16g s 

m— 1 
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The corresponding 6u/fc primary operators are in the form 

ex P ^(Q°$o + Q°e° + q c $ c + Q c e c + q°$° + g°e° + q s $ s + g a e,)/(2\/2)], 

with the scaling dimensions 

(<3c) 2 , ffc(Q°) 2 , (Q c ) 2 , g c (Qc) 2 , (Q°) 2 , 5 S (Q°) 2 , (Q,) 2 , s,(Q,) 2 



165c 



16 



165c 



16 



165s 



16 



165s 



16 



(6.16) 



(6.17) 



(6.18) 



r 



2. boundary states 

Following Cardy, the boundary conditions can be rep- 
resented in terms of the corresponding boundary states 
upon modular transformation. Without going in to the 
detail, we construct the boundary states for the case of 
two quantum wires. 

Because of the conservation of total charge and spin, 
it is natural to impose the Neumann(N) boundary con- 



dition on the center of mass modes, ^c(s)- Since the 
N boundary condition on $° implies Dirichlct boundary 
condition on the dual field 0°, the winding of 0° along 
the boundary should be zero; hence the corresponding 
winding number Q c i s \ = 0. Because of the gluing condi- 
tions, the quantum numbers Q°(s) arc now restricted to 
even numbers and follow some extra constraints. The N 
boundary state for both center of mass modes is 



oo oo oof 

|JV°) = G^exp[^a^ t a^ t ]exp[^a^ t a^ t ] £ expHQ°0° iC ] exp[-iQ° s Q JQ° c , 0) ® 0), (6.19) 

"=1 n=l Q° (s) =-oo 



where prime over the summation indicates the gluing 
conditions, and G% is the ground state degeneracy which 
will be fixed by Cardy 's consistency condition [28l.l30l.[33|. 
The vacuum states are denoted by the winding number 
\Qfti Qfi) °f each independent boson field. For simplicity, 
the phases 0[j c , g % correspond to the applied voltages in 
the wires and will not affect the boundary physics; hence 



we will take these phases as zero. 

Despite the fixed boundary condition for the center of 
mass modes, the other degrees of freedom can have differ- 
ent boundary conditions. In the case of a junction of two 
wires, there are two possible boundary states, Neumann 
and Dirichlct boundary conditions on the ip field (N and 
D boundary states). Similar to the case of the center 
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of mass mode, the N boundary state for the dynamical fields 5? c (s) is given by 
I 



\N c{s) ) = Gjv, c(s )cxp[^a^ c(s) a* c[s) ] ^ |Q C ( S ),0>. 



(6.20) 

The D boundary state implies that the winding of $ c (s) along the boundary is zero. As a result, the quantum number 

Qc(s) — an d Q c (s) is restricted to even integers with more constraints from the gluing conditions. Hence, the 
boundary state is given by 



\D c(s) ) = G DAs) exp[- a n,c(s) a n,c(s) 1 Yj I ' 
n=l Q c ( 3 )=— oo 

where the minus sign inside the exponential function comes from the identification, 4> R = 



(6.21) 



C. 



mensions of primary boundary operators. However, the 
gluing conditions complicate the calculation of the parti- 
In principle, we have everything we need to compute tion function. If the charge and spin arc conserved, there 

is a simpler way to find the dimensions of the boundary 
operators. Since the boundary state of ^(s) * s Neumann, 
(6.22) only the Q^ s ^ = states occur. Hence the constraints of 
Eq. (|6.15[) reduce to 



3. projection of center of mass modes 

ciple, we ha^ 
the partition function 



(B\e- lH v \B) 



with a given boundary state and gluing conditions. Upon 
modular transformation, one can read off the scaling di- 



J 



i- Qc = Q Q S = Qc = Q s = Qc = Q s 

2. Q° C +Q° S =0 

3. Q° C + Q S +Q C = 

4. Q S + Qc = 

5. Q Q C + Q Q S + Qc + Qs + Qc 
From the first constraint, we obtain 

Q° c = 2fi° c , Q° s = 2n° s , Q c = 2n c , C 
Here, all arc integers and follow the constraints below 

2. n° c + n° s = 

3. h° c + n s + h c = 

4. n s + n c = 



Q 

= 2n s , 



= (mod2) 
(mod4) 
(mod4) . 
(mod4) 
(mod8) 



(6.23) 



5. rvi + iri + n c + n s + h c + h s 







Q c = 2h c , Qs = 2h s 



(mod2) 
(mod2) 
(mod2) 
(mod4) 



.24a) 



.24b) 



where the numbering of equations corresponds to that of Eq. (|6.23p . There are four sets of possible parameterizations 
of n M that satisfy all the constraints 



(6.24c) 



'(a) 


n° c 


= 2m°, 




= 2m°, 




n c 


= 2rh c , 




n c 


= 2m c , 




n s 


= 2m s , 




n s 


= 2m s 




(b) 


n° c 


= 2m°, 


n° s 


= 2m°, 




h c 


= 2m c H 


- 1, 


n c 


= 2m c H 


-1, 


h s 


= 2m s 4 


- 1, 


n s 


= 2m s - 


- 1 


(c) 


n° c 


= 2m° + 1, 


fi° s 


= 2m° - 




h c 


= 2m c H 


-1, 


n c 


= 2m c , 




h s 


= 2m s H 


-1, 


n s 


= 2m s 




Ad) 


"S 1 


= 2m° + 1, 


n° s 


= 2m°- 


HI, 


n c 


= 2m c , 




n c 


= 2m c H 


-1, 


h s 


= 2m s , 




n s 


= 2m s - 


h 1 



with the constraint rhP c 
(mod 2) for each set. 



In order to disentangle the Hilbcrt space of ^(s) f rom 
the total Hilbert space, we denote a state with quan- 
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turn numbers Q c ( s ) and Q c ( s ) m the Hilbert space of 
^c(s) by q }. Note that this state is propor- 

tional to \Q c (s),Qc(s)}- Although only \tpQ, ) is non- 
zero in the N state and only |V'o.Q c(s) ) is non-zero in 
the D state case, we would like to consider the general 



case. Using Eq. (|6.24ap and Eq. (|6.24cp and considering 
that m°, s) = 2l c{s) or ^c(s) = ^c(s) + 1j the boundary 
sates corresponding to the center of mass modes can be 
generally written as 



K c{s) ) 



(4 5 c( s) ) 1/4 exp[^ 



0,L t o,R t-i 

a / x a , \ 



oo 

(^( S )) i/4 ex P E/i;a:j;] e i^+^o), 

l c{s) ez 



n=l 
oo 



(4 5c(s )) 1 / 4 ex P [Ea^( s) t <c W t ] E K)+2,0>, 



~n,c(s) "n,i 
»=1 i c(s) 6Z 



(4 5c(s) ) 1 / 4 e Xp [E<i s) t <^ ) t ] E l«Sw + «.0>, 



(6.25a) 
(6.25b) 
(6.25c) 
(6.25d) 



where combinations \Nq and 1-^2 c ( s )) correspond to the sets (a) and (b) in Eq. (|6.24cp and combinations \N® c , s y 
and 1-^3 c (s)) correspond to the sets (c) and (d). Hence, the most general boundary state can be written as 



\B) = (|< c > ( 


s i< s ; 


+ l^2°c> 


8 l^ 2 ° s > 


) ® E (\^m c ,4m c ) 8 


IV^m* ,4m s ) ) 




+ (\Kc) c 


3 l^ 2 ° s ) 


+ l^2,c) « 




® E (I^4m c ,4m c ) ® 


|V'4m,,4m,}) 




+ (l<c) < 


3 K) 


+ l^.c) C 


5 1^2%) 


® E (h/ , 4m c +2,4m c +2) ® |^4m s +2.4m s - 


h2» 


+ (l<c) < 


3 l^ 2 ° s ) 


+ Kc) C 


5 |< s ) 


® E (h/ , 4m c +2,4m c +2) ® V'4m s +2.4m s - 




+ (K,c> c 


3 |JV?,.> 


+ |JVa, c) C 


3 |iv 3 ° s ) 


® E (1^477^+2,4^,) 


<£> fe,+2,4m,» 




+ (l<c> C 


3 |iV 3 ° >s ) 


+ Kc) C 




® E (IV'477 lc +2,4r77 c ) 


<g> |-04rh s +2,4m s )) 




+ (Kc> c 


3 |< s ) 


+ l< c ) C 


3 |iv 3 ° s ) 


(8 E (l^4m ,4m +2) 


® fe„4m s +2» 




+ (l<c) < 


3 |iv 3 ° s ) 


+ 1^3° c ) C 


5 l< s ) 


® E (IV , 47T lc ,4m c +2) 


® |V'4m„4m,+2}) 





(6.26) 



where the prime over the summation indicates the constraint rh c + m c + rh s + m s = (mod 2) for each term and the 
double prime over the summation indicates the constraint m c + m c + fh a + m s = 1 (mod 2). 
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Using Eq. (|6.22[) . the diagonal partition function is given by: 

/ 

^ ' (^4m c ,4m c *Zu±m s Am^) 
ll 

+ (^Ng c Zj^a + Z N o Z N o ^) ^ ^ (^4m ci 4m<, ^4m„4m s ) 
+ (Zjvo ^ Z^o s + 2jvO c Zjyo ^ 
+ {Z N O c Z N s + Z N o^ Z N o^ 
+ ( Z Nl c Z N b + Z N a^ Z N o^ 
+ ( Z N° C Z N o s + Z N O c Z N a s 
+ { Z N° C Z N a s + Z N c Z N o^ 
+ ( Z N° C Z N° B + Z N° C Z N o 



(Z4m c +2,4m c +2 Z^fh,+2,4m B +2) 

ii 

{Zifh c +2.im c +2 Zif ns+ 2Am B +2) 

I 

(Z4m c +2Am c Znh, +2Am, ) 

// 

(Z4m c +2,4in c Z4 7 f ls+ 2,4m s ) 

I 

(Z4m e; 4 mc +2 •Z'4m si 4„i 3+ 2) 
(Z4m e; 4„ lc +2 Z4m si 4„i s+ 2) 



(6.27) 



Here, Z^o ^ ^ are diagonal partition functions in the $°(s) Hilbert space while the ( } q ( } are diagonal partition 
functions in the $ c (s) Hilbert space: 



\e~ lH e lib- ) 
«(«)««(«) "Vcw.Vcw lr Qc (s ),Qc( ii )' 

The partition functions corresponding to the center of mass mode modes are given by 

1 



Z N o 



0,c(») 



«7(?) 



c(«) 



Z N o 



Z N o 



3,c(s) 



C( S ) 



77(g) ^ 9 
1 



(Qc (s )) 2 /(8fl c(s) ) 



y^g2 gc(a) (Z° (3) + l/4) 2 = _J_ y- e -*f<^ (s)9 (Q° (s) ) 2 /(89c( S >) ; 



c(«) 



Vjgcfrj \- ~2 3c , s ,(F°,., + l/2) 2 _ 1 V p -"Qc(a) ff (Qc( s )) 2 /(85c( s) ) 



c(») 



«7(9) 



r^(.)S„,+3/4) 2 = l. y e -^Qc (a)(Z (Qc( S ,) 2 /(89 C(s ,) 



(6.28) 

(6.29a) 
(6.29b) 
(6.29c) 
(6.29d) 



where the modular transformation has been performed 
at the second equality and 

q = e~^, q = e~^. (6.30) 
Here, we introduce the Dedckind ?^-function: 



n 



(6.31) 



which comes from the oscillator modes of the boundary 
state. One can show that each set of independent os- 
cillator modes contribute a factor of l/r)(q). Moreover, 
the modular transformation of the Dedekind 77-function 



is given by 



(6.32) 



The scaling dimensions in Z N o correspond to the 

i,c(«) 

boundary operators that include 3>°(s) an< ^ therefore do 
not conserve the total charge and spin. For instance, 
a nontrivial dimension, (Q®^) 2 /(8<7 C ( S )), corresponds to 

the vertex operator e 4( ^ c ( s >* c < s >^ 2v/ ^ ) . Since we only con- 
sider the perturbations which conserve the charge and 
spin, those vertex operators should not appear. The term 
1/77(5) corresponds to irrelevant or marginal boundary 
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operators which come from the derivatives of an< ^ 
have integer dimensions. 

Thus, we shall concentrate on the boundary operators 
involving only the dynamical fields < & c (s) and replace 



mass modes from the full boundary states. The original 
gluing conditions can be reduced to a set of constraints 
on Q c (s) an d Q c (s) only, which will be refer to as the 
"Reduced Gluing Conditions" (RGC). 



Z N o -> 1. (6.33) 

i,c(s) 

The effective partition function becomes 

ZbB ► 2 S ' (^4m c ,4m c -^4rfi s ,4m s ) 

+ 2 (^4ra c +2.4m c +2 Zim s +2 ,4m s +2) 



(6.34) 



Here we outline the procedures for obtaining the scal- 
ing dimensions of boundary operators from the reduced 
partition function. First, we eliminate the center of mass 
modes from the original boundary state in Eq. Q6.26p to 
obtain the "Reduced Boundary State". Then, we derive 
the nontrivial RGC for Q c i s ) and Q c ( s ) from the original 
gluing conditions. Finally, we compute the reduced parti- 
tion function upon a modular transformation and obtain 
complete spectrum of the scaling dimensions of boundary 
operators for a given boundary state. A similar reduc- 
tion procedure can be applied to the case of junctions of 
three quantum wires for spin- 1/2 electrons. 



m£Z 

+ 2 {Zim c +2Am c ^4ni s +2,4m s ) 

m£ Z 

+ 2 (^4m c .4m c +2 ^4m,,4m,+2) 

mE Z 

Observe that this reduced partition function can be ob- 
tained by simply eliminating the charge and spin center of 

I 

C. Partition functions and scaling dimensions of boundary operators for junction of two quantum wires 

In this subsection, we apply the method described above to compute the reduced partition functions and scaling 
dimensions of boundary operators involving the dynamic field & c {s) for given boundary conditions. 

1. NN Boundary State 



The reduced NN boundary state can be written as 

oo oo / 

\NN) = g NN exp[£ a^ c + ]T a£ a£] ]T \(Q C , 0)) ® \(Q S , 0)>, 



(6.35) 



where gNN is the ground state degeneracy and prime indicates that the RGC should be obeyed. The corresponding 
partition function for the NN boundary state is given by 



J N N ,N N 



= (NN\e~ lH f> \NN) = ip^exph(^(Q c ) 2 + ^(Qs) 2 )} 



8/3 



8/i 



(6.36) 



r 



Now, we shall discuss how to obtain the RGC. Recall- The conditions a and b imply that h 



and fi r have 



ing the gluing conditions in Eq. (|6.15p and using Q° 
Qc(s) = 0, one concludes that Q c and Q s are even. Hence, 
Qc!s) = 2tt. c (s)- Still, we should keep in mind that Q° c and 
Q° s are even and contribute to the reduced gluing condi- 
tions. In terms of n^, the gluing conditions are reduced 
to 



b. h 



n c 



c. hi + h c 



: 

= 

h°„ 



= 



(mod2) 
(mod2) 
(mod4) 



3.37) 



the same parity. In addition, one can conclude that h s 
has the same parity as h c and n®/ s \ due to the condition 
c. Hence 



= 



(mod2) 



(6.38) 



becomes the only gluing condition for the variable h c ( s y 
Naively, one may expect that condition c should provide 
other constraints. However, it can be shown that the 
extra condition is redundant. We justify this statement 
below. 

First, because all the winding numbers, n^, are of the 
same parity, we can set either = 2m ^ or = 2771^+1. 
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For either case, the condition c leads to 

m° c + m c + m° + m s = (mod2) (6.39) 

However, any combinations of integer rh c and m s have 
relative sets of fhP c and fhP s which satisfy this gluing con- 



dition for m p . Moreover, since the partition functions 
corresponding to center of mass modes provide an equal 
constant, set to be 1, for any combinations of J7W S ), this 
condition effectively provides no constraint on h c r s y 

Now, the partition function, Eq. (|6.36p , becomes 



z NN , NN = fpE ex p[-(^^) 2 + lj $&>n . ( 6 - 4 °) 

in terms of variable n M with the reduced gluing condition, Eq. (|6.38[) . Further, the RGC leads to two separate sums 

Z NNNN = 4^{ T tfWJ'+W*'?) + Y ^(2^+i) 2 +¥(2™ s +D 2 ) } (6 .4i) 
,yl ' m c(s) ez m <a) £Z 

without any constraint on rh c and rh s . Upon the modular transformation, the partition function becomes 

Z NN , NN { q ) = 4=^2 E ^ + ^ X [l + e -Mm c+ms)] (6 42a) 

= £ q ^^f . (6.42b) 

Hence, the spectrum of boundary operators can be read off and the leading order perturbations have the scaling 
dimensions 

112 2 , 

^ + ^, -, (6-43) 

2.9c 2g s g c g s 



which agree with the previous results obtained using the DEBC method in Sec. IIVI Moreover, the ground state 
degeneracy, useful for determining the stability of the fixed point, can be found by Cardy's consistency condition as 
9nn = (ff c .9s/4) 1/4 . 



2. DD Boundary State 

The reduced DD boundary state can be constructed as 

oo OO / 

\DD) = g DD exp[- E <}c<% - E a n,s<}s\ E K°> <W> ® K ' &)>> ( 6 ' 44 ) 

n— 1 n— 1 

using Eq. (|6.2ip with ground state degeneracy gDD- Again, the prime over the summation implies the reduced gluing 
conditions. The partition function corresponding to the DD boundary state is given by 

Zdd.dd = (DD\e- l ^\DD) = £ C xph(^-(Q C ) 2 + ^(Qs) 2 )} , (6.45) 



Since Q®r s \ = Q c (s) = 0, the condition 1 of Eq. (|6.15|) numbers become 
) = 2 ^c( s ) and Qc(s) = 2n c(s) for integer 

=> crliiitiff rnnrlitinns in t.prms nf t hp niiflnt.iim 

b. h° c + n s = (mod 2) 



leads to Q° ( ) = 2n° ( s and Q c i s \ = 2?i c(s) for integer _ 

c . w . . cw ... V r , a. + rvi = (mod 2) 

71^. Hence the gluing conditions m terms of the quantum 



(6.46) 

c. n c + n s = (mod 2) 

d. 7i" + n c + n° + n s = (mod 4) 
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Here conditions a-c imply that n c ( s ) and "°(s) have the with the constraint rh c + rh s = (mod 2). Therefore, 

same parity. Condition d provides no constraint on the the scaling dimensions of the leading order perturbations 

quantum number n c ( s ) by a similar reason as in the case can be read off as follows 
of the NN boundary state. 

Upon the modular transformation, the partition func- — — — , 2g c , 2g s . (6.48) 
tion corresponding to the DD boundary state can be eval- 2 

U.atcd as 

This is in agreement with the previous results obtained 

2 / ' using the DEBC method in Sec. IIVI Finally, the ground 

Z D D.DD{q) = 9 2 dd^^Y. <Z* ( " lc)2+ * ( ™ s)2 (6.47) state degeneracy is given by g DD = (l/4g c g s )^ 4 . 

I 



3. ND Boundary State 

We first construct the reduced boundary state for the ND boundary condition 

oo oo / 

\ND) = g ND exp[]T a^< f c - E a nWU E °)) ® K°> 0*))' ( 6 ' 49 ) 

71 — 1 11— 1 

where the qnd is the ground state degeneracy. Now, the corresponding partition function can be computed as 

Znd.nd = (ND\e- lH ? \ND) = ^ £ e xp[-(^(Qc) 2 + ^(Qsf)] , (6.50) 

with proper reduced gluing conditions indicated by the prime over the summation. 

For obtaining the reduced gluing conditions on Q c and Q s , we first observe that Q^ s y Q c and Q s are even integers. 

Hence, Q°( s ) = 2n^ s y Q c = 2h c and Q s = 2n s . After some algebra, one finds the following reduced gluing condition, 

h c = n s = (mod 2) (6.51) 
Upon the modular transformation, the partition function becomes 



Z ND<ND { q ) = yi^L^^C-) 2 ^^) 2 , (6.52) 

without any constraint on m c and rh s . Thus the scaling dimensions of boundary operators for the leading order 
perturbations are given by 

f , ^ ■ (6-53) 
2 2g c 



Again, this matches the results from the DEBC method in Sec. IIVI The ground state degeneracy becomes {g c /9s) 



1/4 



4. DN Boundary State 

First, we construct the reduced boundary state for the DN boundary condition 

00 00 / 

\DN) = g DN cxp[- J2 < f c <c + E E K°> Qc)) ® \(Qs,0)) (6.54) 

n— 1 n— 1 

with gDN defined as the ground state degeneracy. Then, the corresponding partition function can be written as 

Z DN , DN = (DN\e~ lH v \DN) = £ eX p[-(^(Q s ) 2 + ^(Qc) 2 )], (6.55) 



I 

with the proper gluing conditions. Again, all relevant 2nP^ s y Q c = 2n c and Q s = 2n s . One can show that the 
momentum quantum number are even. Thus, Q c r s \ = 



2G 



only reduced gluing condition for n c and h s is 



= h s = 



(mod 2). 



(6.56) 



Upon the modular transformation, the partition function 
becomes 



J DN,DN 



(<?) 




(6.57) 



with arbitrary integers, rh c and m s . Hence, the scaling 
dimensions of boundary operators for the lowest order 
perturbations can be read off as 



9c 
2 ' 



1 

2g~ s 



(6.58) 



This is in agreement with the previous results from the 
DEBC method in Sec. HVl The ground state degeneracy 
is given by (5 s /g c ) 1/4 - 

In summary, we have computed explicitly in this sub- 
section the spectrum of boundary operators and the 
ground state degeneracy for a junction of two quantum 
wires. The use from the onset of the U c (l) x U s (l) sym- 
metry, corresponding to the total charge and spin con- 
servation, leads to the reduced boundary states and re- 
duced gluing conditions; these reduced relations largely 
simplify the computations. It is worthwhile to emphasize 
that this simplified scheme can be applied only to a sys- 
tem with the conservation of total charge and spin and 
without resonant tunneling. 



VII. JUNCTION OF THREE QUANTUM 
WIRES FOR SPIN-1/2 ELECTRONS-BCFT 

In this section, we will apply the technique developed 
in the previous section to the case of a junction of three 
quantum wires for spin-1/2 electrons. We shall first 
derive the gluing conditions for a convenient basis and 
project out the center of mass modes of the charge and 
spin degrees of freedom. Then, the reduced partition 
functions for given boundary conditions and the scaling 
dimensions of boundary operators will be computed and 
used to determine the stability of the fixed points. 



and a corresponding set for 9 fields. 

The momentum quantum numbers of the total charge 
and spin modes follows the relations 



Q°c(s) 

Qc(s) 



= Ql 



c(s) 



Q 2 



c(s) 



Qc(s)- 



(7.2a) 
(7.2b) 



Again, due to the conservation of total charge and spin, 
the N boundary condition should be imposed on the $°(s) 
field. The corresponding boundary state always has the 
quantum number Q®r s \ = 0. Hence, it is convenient to 
parameterize the vector of integers, Q % c i s y as: 



(Ql(s)>Qc(s)>Qc(a)) 

ml {s) (0,l,-l) + ml (s) (-1,0,1), 



(7.3) 



where to 1 ! 2 ? G Z. Consequently, to 1 



~Ql(«\- O n tne other hand, Q 



Q c (s) 



and 
\ can De 

nonzero. So it is convenient to parameterize the Q' c ^ 
as 



o 

c(s) 



(Ql 



c(s) 



EE n° 



(s) (l, 1, 1) - n l c(s) (0, 1, 1) - 7i 2 (s) (l,0, 1). (7.4) 



Therefore, ril,,^ can be expressed in terms of Q l c , s \ as 



1 

c(s) 



fc(s) 



fc(s) 



-Q 3 



c(s) 



(7.5) 



Hence the "*(s) P rov ide a representation of quantum 
numbers in an alternative basis. 

Now, let us investigate the corresponding gluing con- 
ditions of the new variables m l , > and n l , From the 

c(s) c{s) 



gluing conditions in Eq. (|6.10[) and the definition of Q° 
and Q®/ s \, one concludes following gluing condition, 



c(s) 



Ql = Ql = Ql = 







(mod 2), 
(mod 4). 



(7.6a) 
(7.6b) 



Because Q®^ = 0, the following gluing conditions hold 



A. Reduced gluing conditions 

Here, we start from the mode expansion of the bosons 
for a single quantum wire, Eq. (|6.12a[) and Eq. (|6.12bD 
with the gluing conditions in Eq. (|6.10[) . and generalize 
the mode expansions and gluing conditions to another 
orthogonal basis: 



c(a) 



1 

71 
1 

71 
i 

71 



(Vc( 8 ) + <P%) + Vc( S )) 



In terms of : 



Ql = Ql = (mod 2), (7.6c) 
Ql + Ql = (mod 4). (7.6d) 

'/ s, the first condition Eq. (|7.6cp becomes 



Q°c(a) = 3 «c(.s) " 2 ("c(s) + = (mod 2). (7.7) 



This implies that = 2Pc( s ) are even integers, where 
p'c(s) are arbitrary integers. By using the gluing condi- 
tions for each of the quantum wires Eq. (|6.10p and the 
(7.1) fact that n^,-, are even, one can prove the following glu- 



i t \ 

c(s) 

ing condition 



m\ (mod 2), 



(7.8) 
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from the relations between TO*^(n*/ v) and Q^^iQl^ 
Combining the gluing conditions Eq. (|7.6dj) and Eq. {77 
a new gluing condition emerges 



Finally, two nontrivial gluing conditions for m l c t s ) an d 
n l , > arise 

c(s) 



P°+Z>° = (mod 2). 



(7.9) 



Q 



1(2) 



Q 



1(2) 



)1(2) 



Ql [2) = (mod 4) 



l ?W+7^( 1 ) +m 2(l) + , 



= 



= 



(mod 4), (7.10a) 
(mod 4), (7.10b) 



where the third equality holds because n° + rP s = 2(p" + p") = (mod 4). 
In summary, we found the following gluing conditions for the n l c ^ and m 



a. n , 

c(s) 

b. ni = ni 



2rf 



c(s)> 

mi 



d. ml + ml + nl + 



mi 

nl = 
nl = 



"c( S ) 

(mod 2) 
(mod 2) 
(mod 4) 
(mod 4) 



(7.11) 



In particular, ^c(s) s disentangle from the rest of quantum numbers and do not have nontrivial gluing condition. 



B. Mode Expansions and center of mass mode projection 

We first define a vector field and a conjugate vector field, representing the dynamical boson fields, 



c(a) 



f$ 1 , , $ 2 , ,) 

V c(s)i ^c(s)/i 



e 



c(a) 



(7.12) 



where j are defined in Eq. ()7.ip . Using the definition of the dynamical field in Eq. (|7.1j) , the periodicity along the 
spatial direction follows 



A6 C(S ) = e c(s) (/3,t)-6 c(s) (0,t) = V2ir(m 1 cis) K 1 +m 2 cis) K 2 ), 
where n 1 ,-^ and m% c ( s ) a re defined in the previous subsection and -Ki(2) and i?i(2) are defined as 

Ki = «i - 4(^1 x 2) = (+1,+^), 

#2 = (-5,-^), 4 = ^=(*2 X *) = (-1,+^)- 

Hence the mode expansion of the two-component boson field becomes 



(7.13a) 
(7.13b) 



(7.14) 
(7.15) 



<I> 



c(s) 



1 00 1 

= V — 



where the integers ^(s) an< ^ m r/il are restricted by the gluing conditions. The corresponding energy is 



c(s) 



^iE^i 2 + fiE4' 2 + i'E^i 2 + |iE4' 



(7.16) 



(7.17) 
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where • • • represents the energy of the oscillator modes. 

Recalling the discussion of the last section for the case of the two quantum wires, the center of mass modes will 
only contribute a constant and can be projected out with a corresponding reduced gluing condition. Observe that 



Q°c(s) = 2(3P° 



c( 8 ) 



Ms) 



) implies that there are three classes of quantum numbers categorized by Q°( s \/2 



(mod 3), we shall introduce auxiliary quantum numbers fc c ( s ) = n c , a ^. +^(«)- 



c(s)l 

Thus a general boundary 



state which is N with respect to $°, i.e = ^' takes the form 



\b)= £ 

fc c(s) = -l,0,l 



E 



(\i>{n\y c ,m l c ,ml)) ® \il>{n l s ,n 2 s ,m l s ,m 2 s ))) (7.18) 



-fc c ( s ) (mod 3) 



where the prime over the summation indicates the gluing conditions. With the gluing condition, the boundary states 
corresponding to the <I> fields are given 



l<) ® l<) = (9 5 c. 9s ) 1/4 £ |2(3p° + fc c ), 0) ® |2(3p° + fc s ), 0). 

p0+ p o=0 (mod 2) 

Hence the corresponding "partial Neumann" partition functions are readily calculated 

V 9cg S \^ ~[g c (3 Pc +k c ) 2 /12+g s (3 Ps +k B ) 2 /12] 



Z N? N? - 



r)(qf 



E 



pO+ p o =0 ( mo( j 2) 



(Qc) , (««) 



(7.19) 



(7.20) 
(7.21) 



with the constraint Q® c + = (mod 2). As in the case of two quantum wires, this part of the partition function 
corresponds to boundary operators changing the total charge and spin of the system. Hence, Z N o N o can be projected 

fcc fcs 

out and replaced by unity. Now, the dimensions of all primary boundary operators involving only the dynamical fields 
can be obtained by using a reduced boundary state, which lives in the reduced Hilbert space of the 2-componcnt 
boson fields, $ c ( s ). Thus, the reduced boundary state becomes 



\B) 



EW(ni.»c.»»i.m?)) ® \^(nl,n 2 a ,ml,m 2 s )}), 



(7.22) 



with the constraints given in Eq. (|7TTTj) . Note that there is no gluing conditions between _P°(s) anc ^ the res ^ °^ quantum 
numbers, hence we can simply ignore ?Os) an< ^ take the rest of gluing conditions as the constraint for the reduced 
boundary state. 



r 



C. Reduced partition functions and the dimensions 
of boundary operators 

We first recall the reduced gluing conditions for the 
quantum numbers ^(s) an< ^ 



c(s) 



trivial scaling dimensions of the operators. We shall in- 
vestigate below the scaling dimensions of boundary op- 
erators for given boundary conditions. 

1. NN boundary condition 



a. 
b. 



m{ 

ni = 



(mod 2) 
(mod 4) 



(7.23) 



where j = 1,2. In obtaining the reduced partition func- 
tions, these RGCs play crucial roles and lead to the non- 



As discussed in the last section, the corresponding 

1(2) 

boundary state of N BC has quantum numbers Q c r s ( — 0- 

1(2) 

This implies m c{ s ) = and the NN boundary state is 
given by 



J 



\NN) 



OO 



_^ a n}s a ri]s\ ^ j 



|(t&t&0 J 0))®|(nJ,T&0 ) 0)), 



(7.24) 
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where the prime over the summation indicates the gluing conditions and gNN is the ground state degeneracy. The 
corresponding partition function is calculated as 



Z NN = (NN\e- lH ?\NN) = ^™ £exp( 



813 



(7.25) 



Because rn 1 ^} = 0, n^^s are even using the condition a of Eq. (|7.23p , hence we conclude that ?\fs) = ^c(s) wn -h 



1(2) 



!(2) _ o^!(2) 



l> I \ 

c(s) 

1(2) . . . 1(2) 

arbitrary integers for h c ^ s y Rewriting the condition b of RGC in terms of h c ^, one obtains 



c(s) c(s) 



h{ + hl = (mod 2), for j = 1, 2. 



(7.26) 



It is evident that h{ and h J c have the same parity. Hence, there are four plausible combinations of [h\,h^,h%h? s ) 
categorized by even or odd integers and listed below 



(e,e,e,e) I 
(o, o,o,o) ... II 
(e, e, o, o) . . . Ill 
(o, o, e, e) . . . IV 



(7.27) 



where e stands for even while o stands for odd. The total partition function is the sum of the partition functions 
corresponding to each combination above. 

In case I, we can introduce new integer quantum numbers ^(s) = ^c(s) ano - compute the corresponding partition 
function 



c(s) 



where the second equality is obtained by multi-dimensional modular transformation. Appendix [B] and 4^ » are 

arbitrary integers. Similarly, one can execute the calculation for case II by introducing h ,-, = ^c(s) Upon the 
modular transformation, the partition function is given by 



(7.29) 



with arbitrary integers i^(s)" ^ or cases m an< ^ IV, the similar parameterization of h J r (^ leads to the following portions 
of partition function 



c( S ) 



ylll 



7 IV 



3 ffLv VqC^IS^gjl'+ah-ISj^l') x p-"(*o+* a 1 ) 



_ 3 g^Af ^^lEi^-r+^iEiti^i 2 ) 



x e 



(7.30) 
(7.31) 



<=(«) 



By adding contributions from each part, the full par- 
tition function is given by 

^jvn - rrr 7.9 29c 293 ■ ('-32j 

Observe that the unit vectors Kj form a triangular lat- 
tice. In general, the scaling dimensions can be calculated 



by finding the length square of linear combinations of Kj 
with the constraints t J c + ti = (mod 2) for both j = 1,2. 
The dimensions of the boundary operators corresponding 
to the leading order perturbation is given 

11 2 2 

A NN : + -, -, (7.33) 

2g c 2g s g c g s 
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which are the same as that for junction of two quantum DEBC met hod. M oreover, the ground state degeneracy 
wires with NN BC and agree with the results from the is g^N = y <?e<?s/3. 



J 



2. DD boundary condition 

1(2) 

The corresponding boundary state with DD BC takes the form with nJJ = 

oo oo / 

\DD) =.9DDCxp[-(X;4! c a5 + E^OEl( ' ' m ^ m ')) ® l(0,0,mi,m;)), (7.34) 

ri— 1 n— 1 

with the proper gluing conditions indicated by the prime over the summation. The diagonal partition function is 

Z DD = (DD\e~ lH v \DD) = ^ £ex P [-^(l| £m^| 2 + U 5>l#,| a )]. (7-35) 

Let us investigate the reduced gluing conditions here. Again, the condition a of Eq. (|7.23p and n ! = lead to 



C ( S ) 

even m c(s)' ncncc m c( s ) = ^ w c(s) ^ or arbitrary integers ^(s)' Then condition b of RGC can be written in terms of 

<4 + wi=0 (mod 2), for j = 1,2. (7.36) 

Similarly to the case of the NN BC, there are four possible combinations of (w^, , uj 2 , ui 2 s ) categorized by even or 
odd integers and listed below 



(e, 


e, 


e, 


e). 


.1 


(o. 


o. 


0. 


o). 


.11 


(e, 


e, 


o. 


o). 


..III 


(o. 


o. 


e, 


e). 


..IV 



k,w s ,u c ,w s ) = < . (7.37) 

(e,e,o,o) III 

[ (o,o,e,e) IV 

Again, we will investigate the contributions to the partition functions of each combination. In case I, we can 
introduce new integer quantum numbers v J ,-. = 2i^ s j and compute the corresponding part of the partition function 

for arbitrary integers ?( s )' Since Ri forms a triangular lattice with the lattice spacing = it is convenient to 

introduce scaled vectors R'i = ^Ri with unit length, \R'%\ = 1. Hence the partition functions in terms of these unit 
vectors are given by 

z i = Qc&fpD Y w\ £^R',| 2 +^| fj^l 2 ). (7 .3 9) 
DD 12 r)(qY ^ H y ' 



For other cases, we can introduce new quantum numbers l° c ^ such that oj 3 , — 2i^,N + 1 for odd integers and 
w c(s) = ^c(s) ^ or even integers. Hence the partition functions for each case are given by 

z u y c y s y DD (m Tl, gjg < I y., Hr>^) -frfe+g+g+g) f7 40a) 





.9dd 


12 


^) 4 


9c9s 


9 2 DD 


12 






9 2 DD 



(7.40b) 



z iv = ycys y DD y> (»iy;,gjp < i 3 +»iy; J gjp < i a ) c -fa(e+g) / 740c ) 

12 77(g) 4 ^ 

*c( S ) 



I 

Now, let us add all contributions and obtain the full with the constraints t° c + t{ =0 (mod 2) for j = 1,2. 
partition function The scaling dimensions of boundary operators for lead- 

Z DD ^iff^± q ^^^ + ^^^ (7 .41) 



31 



ing order perturbations can be read off directly from the 
partition function 



A 



DD 



9c 
6 



G 



25c 
3 



2?s 
3 



(7.42) 



Again, the results agree with that obtained by the DEBC 
method. Finally, the ground st ate deg eneracy of the DD 
boundary state reads Qdd = \/3/ 1 g c 9s- 



can be calculated by exchanging the charge and spin sec- 
tors of the ND case.) The boundary state corresponding 



to the ND BC can be constructed with m 



1(2) 



and 



1(2) 







3. ND and DN boundary condition 

We will focus on the ND BC first. (The scaling dimen- 
sions and partition function corresponding to the DN BC 



\ND) = (jND cxp[+ £ i£t cS flt _ £ a£t aS flt ] £ n l o, 0)) ® |(0, 0, m s \ m*)}, (7.43) 

n— 1 n— 1 

with proper constraints and ground state degeneracy gND- Since the gluing condition a in Eq. (|7.23|) implies that 

ni and ml' 2 ' arc even when m^ 2 ' = ns' 2 ' = 0, we shall parameterize = 2hc and m^ 2 ' = 2ojs ( ' 2 '. The 
partition function thus can be computed, 

Z ND = (ND\e~ lH f> \ND) = ^ £ exp[-^(|| £ KM 2 + j\ £ «>%\% ( 7 - 44 ) 

I 



with the gluing conditions K{ + u) J 8 = (mod 2) for 7 = 
1,2. 

After some algebra, similar to the case of the DD and 
NN boundary conditions, one obtains the full partition 
function 



Z 



ND 



9s9 ND ( — 



Ej tjKjl , gal E, *lfl',l 2 



(7.45) 



with the constraints t J c +t 3 3 = (mod 2) for j = 1, 2. Then 
the leading order boundary operators have the scaling 
dimensions 

a 1 , * 2 2 9s ,„ A( ,, 

Aatd : ^— + -7-, — , -5-, (7.46) 
2.9c 6 g c 3 

which match the results from the DEBC scheme. More- 
over, the ground state degeneracy reads g^D — \/ 9cj 9s- 
Due to the symmetric structure between the charge 
and spin part of the boundary states, one can exchange 
c <-> s in Eq. (|7.45[) and the corresponding constraints to 
obtain the partition function given the DN BC. Then, 
the dimensions of boundary operators of leading order 
perturbations with the DN BC are given by 



A 



DN 



(7.47) 



2.9s 6 ' g s ' 3 

The results agree with that of the DEBC scheme . Fur - 
ther, the ground state degeneracy reads qdn = \J 9s 1 9c- 



4- X±X± and x±Xt boundary conditions 

We will first study the case of X+X+ boundary condi- 
tion. The X-X- boundary condition can be calculated 
in the same manner. Indeed, one can show that the rele- 
vant scaling dimensions have exactly the same structure 
for both x±X± boundary conditions. Then, we will com- 
ment on the case of X±Xt- 

Generically, a comformally invariant boundary condi- 
tion can be expressed as 



TR r, TL -.i r> I cos £ — sin£ 
b" = R^, withi? c = . ; 

\ sm£ cose; 



(7.48) 



Since the winding along the boundary can be written as 



AG, 



A0, 



^(V9cjs^^c(s 



y/9c(s) 

A6 c(s) 



2—' 



, (7.49a) 



, (7.49b) 



arbitrary conformal boundary conditions satisfy 

-. A@ c ( s ) AO c f s ) 

y/9cJ^A^ cis) 7 =±± = i?c( v ^yA$ c(s) + nr ±i ). 



\/9c(s) 



\/9c~{s) 
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Observe that AO 



c( S ) 



for the NN BC leads to £ = 



and A$ 



c(«) 



for the DD BC leads to £ = tt. 



For 
and 



£ 7^ 0, 7T, it is clear that we need both A<5 C ( S ) 
AO c ( s ) nonvanishing. Moreover, is a rotation matrix 
which preserves the length of the vector, hence A$ c ( s ) 

and A© c ( s ) have to be mutually orthogonal. For satis- 
fying this constraint, the general oscillator vacua can be 
constructed from 

\(a c T}l,a c riZ,p c ril,0 c r£)) <g> \(a s rjl,a s r] 2 8 ,0 s T]l,0 s r] 2 s )), 

1(2) 

for arbitrary integers V c ls) w ^th proper gluing conditions 



J 



on a and j3. For instance, the condition a in Eq. (|7.23|) 
leads to a c ( s ) = @ c (s) (mod 2). Moreover, the chiral ro- 
tation angle is fixed by (a c ( s ), /3 C («)) as 



tan 



£ c (s) _ /5 c (s) v3 



(7.50) 



Notice that the chiral rotation angle is "quantized" . 

As shown in Ref . 0] , the x+ fixed point corresponds to 
the choice of a c ( s ) = f3 c (s) — 1 while the x— fixed point 
corresponds to the choice of a c ( s ) = — /3 c (s) = 1- Hence, 
the boundary state corresponding to the X+X+ boundary 
condition can be constructed as 



\X+X+) = 9 X 



+x+ 



exp[]T <t c . R( jM c + J2 <t . ^«£tj £ r^ 2 , r,l V 2 )} ® {(rjlvlvlVs)}, (7-51) 



71=1 



ra=l 



where the prime over the summation indicates the gluing conditions. The diagonal partition function is given by 



?/(?T 



/3 4 
tt/ 



(7.52) 



r 



where the second equality hold because Ki ■ Rj = ey and 
the new vectors arc defined as 



A l(s) = V9^) 



Ri 



K, 



2V2 y^7 { 



(7.53) 



These new vectors form a triangular lattice with lattice 
spacing 



14 



/ 9c(s) 



1 



=«' y e 2<? c(s) ' 

Using Eq. (|7.23p . there is only one constraint 
vi = Vi (mod 2), for j = 1,2. 



(7.54) 



(7.55) 



Therefore, there are four possible combinations of the 
winding numbers classified by even or odd integers which 
are listed below 



{vl,vl,v 2 c ,v 2 s 



(e, 


e, 


e, 


e). 


.1 


(0, 


0, 


0, 


o). 


.11 


(e, 


e, 


0, 


o). 


III 


(0, 


0. 


e, 


e). 


IV 



(7.56) 



Again, parameterizing the even and odd variables as 
i] = 21 and 77 = 21 + 1, respectively, we are able to 
compute the diagonal partition function corresponding 
to each set. Upon a modular transformation, the partial 
partition function for each case is given by 
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Z x+x+ 



7ll 



X+X+ 



7 IV 

X+X+ 



(3- 


h 5c 2 )(3 4 


<? 2 ) ??(<z) 4 




12 g c g s 


<? 2 

■^X+X-f- 


(3- 


K5 C 2 )(3 4 


-<? s 2 ) ^) 4 






. 2 

^x+x+ 


(3- 


h 5c 2 )(3 4 


■5 2 ) ??(<z) 4 






9 2 


(3- 


h 5c 2 )(3 4 


-5 2 ) v(q) 4 



q \ E 3 * 3 X j | 2 +| Ej *i^l 2 g-Mtl+tl+tl+t 2 ,) _ 



y g iE;,^r+iSi*^i 

E 

< S ) e 

V g l Ej * J cM>i| 2 + l E, *i W-f 1 3 e — i-rr C*o +*2 ) , 

E 



g i e 3 tj^'r+i Ej * j s M > s 3 i 2 e -^(*c+*') 



(7.57a) 
(7.57b) 
(7.57c) 
(7.57d) 
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where W^' and Wl are the dual vectors of A 3 C and 
respectively, and form two sets of dual triangular lat- 
tices with different lattice spacing, J '^^3 and ^^2+3 j 
respectively. With the identification of the ground state 



degeneracy as g x+x+ 
function is 



4Sg c g s ■ 



the full partition 



1 



J x+x+ 



iE< 



JE.^'I'+IE, tiwi\ 2 



with the constraints P c + P s = (mod 2) for j = 1,2. 
Hence the dimensions of the boundary operators can be 
read off directly and the leading order perturbations have 
the dimension 



A 



25s 



X+X+ • 2 
i/c 



(7.59) 



The only difference of the boundary state of the X-X- 
boundary condition comes from the parameterization of 
a c(s) = —Pc(s) = 0. Consequently, the partition func- 
tion is similar to the case of the X+X+ BC with a minor 
variation, 



c(s) 



.4' 



c(s) - y/9c(e) 



R, 



2V2 y/2g c ( s ) 



(7.60) 



in Eq. (|7.52[) . Also, the A' c ^s form a triangular lattice 
with lattice spacing equal to \A 3 , J. Moreover, the con- 



straint, 7)1 



(mod 2), still holds. Hence, we 



can conclude that the partition function upon the mod- 
ular transformation has a similar structure to the X+X+ 
case. Hence the leading order perturbations have the 
same scaling dimension. 

Since the scaling dimensions of the boundary opera- 
tors are the same for the cases of x±X± BCs, one may 
wonder if the x±Xt boundary conditions will have simi- 
lar behaviors. However, one can see from the argument 



below that the scaling dimensions of the boundary oper- 
ators are rather different. We have learned that the non- 
trivial scaling dimensions come from the gluing condi- 
tions. Even in the case of the same dual lattice structure 
and spacing, the gluing conditions may provide nontriv- 
ial constraints. In the case of the x±Xt BCs, the only 
gluing condition ff c ! 8 ^ = (mod 4) will not lead to any 
constraint on the integer As a result, the leading 



(7.58) 

order perturbations have scaling dimensions 



2ff c 



2ff s 



j X±Xt ■ „2 



(7.61) 



which is always smaller than one for any g c ( s ) and leads 
to an instability of the fixed point. So, we conclude that 
the x±Xt fi xe d points are not stable. 



5. DaDa boundary condition 

From the results of the DEBC approach, the asym- 
metric boundary conditions could be stable in some re- 
gions of the interaction parameter space. Hence, it be- 
comes important to construct the corresponding bound- 
ary state for checking the instability of the DaDa bound- 
ary condition. Without loss of generality, we choose to 
impose Dirichlet boundary condition on the dynamical 
field, s = (<ph s -\ - </? 2 (s) )/V2, between first and sec- 
ond wire, and Neumann boundary condition at the third 
wire. Indeed, this set of boundary conditions is equiva- 
lent to having D BC on Qh, and N BC on $%r s y Using 
the parameters defined in this section, the D BC of the 



dynamical field leads to n* 
BC of the VcU) l ea ds to mVj 



1 



c(s) 



Ms) 



while the N 



1 i \ 

c(s) 



*c(s)- 



Hence, 



the reduced boundary state can be constructed as 
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\D A D A ) = g DA D A exp[J^ a^ c ■ RaS^ c + E a^ s ■ RaS^ s ] E \(n c , n c , m c , m c )) ® |(n s , n s , m s ,m s )), (7.62) 

n— 1 n— 1 

with the rotation matrix and the appropriate gluing conditions written as 

Ra -- 



-1 
1 



a. n c 

b. n r 



lis 

lis 



-- m s 
m. 







(mod 2) 
(mod 4) 



The partition function is given by 
Zd a ,d a = {D A D A \e- lH e\D A D A ) 



4 1 2 



2ttZ , 1 



£^ 2 )]exp[-^ (l i-K£i? ; 



'^4 Z^ L *^ W 



12 n P[ (3 y Ag s 



r 



4 1 2 



(7.63) 
(7.64) 



12 



where we use |£^l 2 = 1 and |£i^i/ 2 | 2 = 1/3 for 
the second equality and the prime over the summation 
indicates the constraints of the integers. 

Now, we can separate n c ( s ) and to c ( s ) into two inde- 
pendent sets, all even or all odd, 



J (a) n c ( s ) = 2/i c ( s ), m c ( s ) = 2w c(s) 
\ (fr) n c (s) = 2/i c ( s ) + 1, m c(s) = 2u> c ( s ) - 

with a constraint 

h c + lu c + h s + lu s = (mod 2). 



(7.65) 



(7.66) 



Hence the partition function can be decomposed in terms 
of new variables as 



J D A ,D A 



< 6 2 9s 6 



(7.67) 



r?(9) 4 

(2^ c + l) 2 , 8c (2h e + l) 2 (2^ s + l) 2 , g s (2fe s + l) 2 , 
q 8g c ' 24 8gs -T 24 I 



Upon the modular transformation, one obtains 



J D A ,D A 



?7(<z) 4 



8 T 8 9s 



with the constraint 

a. t c = t s =t c = t s (mod 2) 

b. t c + t s + t c + t s = (mod 4) 



(7.68) 
(7.69) 



Then the scaling dimensions of the boundary operators 
corresponding to the leading order perturbations can be 
read off, 



A 



D A D A 



9s 



T + ~2~' 2 {9c 



9c 3 
2 2 5s ' 
2ff c , 2g s 



2 



3 

2g~c 



(7.70) 



These results agree with the conclusions from the DEBC 
method. Moreover, the ground state degeneracy is 
9d a d a = 1/V3. 

We conclude this section with the discussion about how 
to determine the instability of the phases by using the 
ground state degeneracy of the boundary states. Recall 
the ground state degeneracy of the different boundary 
states: 



9nn = 
9nd = 

9x+x+ 



9c9s 



3 ' 



9DD 



9c9s 



9DN 



(7.71) 



(3 + g 2 )(3 + g 2 ) 
48g c g s 



9D A D„ 



= 1/V3. 



The universal non-integer ground state degeneracy 9b c b s 
always decreases under the renormalization from a less 
stable to a more stable fixed point in the same bulk uni- 
versality class HH (g-theorem). If there are two stable 
phases for a given value of interaction parametres (g c , g s ), 
then an unstable fixed point must lie in between these two 
stable points (with a value of gB a B, that is larger than 
those at the two stable fixed points). While one cannot 
resolve to which fixed point one flows in these overlap- 
ping regions of stability (it dependes on the strengths of 
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FIG. 9: [Color online] The areas with red, green, blue, or- 
ange, purple and white indicate the regions in parameter 
space where the ground state degeneracy corresponding to 
the NN, DD, ND, DN, x±X± and D A D A boundary condi- 
tions, respectively, is the minimum. The borders with black 
lines are determined by the condition that the ground state 
degeneracies for adjacent regions corresponding to two differ- 
ent boundary conditions are equal. 

the bare couplings), it is instructive to look at the phase 
boundaries as determined from the condition of minimal 
boundary entropy. The phases with minimum ground 
state degeneracy gB c B s [computed from Eq. (|7.7ip ] for- 
given (<? c ,<?s) are shown in Fig. [9] Comparing with the 
phase diagram proposed based on the scaling dimensions 
of the leading order perturbations in Fig. [JJ the border- 
lines in Fig. [5] are located where there is overlap between 
two or more phases. 

VIII. CONCLUSION 

In this paper we studied the problem of a junction of 
three quantum wires for spin-1/2 electrons connected by 
a ring, through which a magnetic flux can be applied. 
The bulk of the wires was formulated as Tomonaga- 
Luttinger liquids with interaction parameters g c and g s 
(for charge and spin sectors, respectively). The problem 
was studied using two different methods: delayed eval- 
uation of boundary conditions and boundary conformal 
field theory. These methods bypass the difficulty that 
normally occurs with the inclusion of Klein factors to en- 
sure the proper fermionic statistics for different species of 
fermions. We reached consistent results for the stability 
of the phases obtained from the two different methods. 

We computed the low energy and low temperature 
charge and spin conductance tensors corresponding to 
the fixed points for the junction as a function of the in- 
teraction parameters g c ( s y These conductance tensors, 
G(g c , <7s), presented in SecfVlfand summarized in Sec HI)) . 



characterize the response of the junction to the externally 
applied voltages. We have presented a simple one-to-one 
relation between the conductance tensor and the rotation 
matrix 1Z that encodes the types of boundary conditions 
in the DEBC method. 

When the Y-junctions are attached to Fermi liquid 
leads, the fixed points are still controlled by the g c and g s 
in the wires. However, the conductance tensor is altered 
due to the contact resistances at the lead/wire interfaces. 
Similarly to what was found in Ref Q, the conductance 
tensor in the presence of the leads is the one determined 
by the appropriate BC (which is controlled by the g c and 
g s ), but evaluated at g c = g s = 1 instead. For instance, 
in the case when the chiral fixed point is the stable one, 
one plugs g c = g s = 1 into Eq. (|2.13p . Interestingly, for 
the chiral fixed point, the switching of the current in the 
presence of the leads is then perfect, circulating the cur- 
rent from one lead completely into one of the other two 
leads. Indeed, this renormalization of the conductance 
tensor is the 3- wire analog of what was found in Tarucha 
et al.'s experiments on Luttinger liquids coupled to reser- 
voirs (34[ and explained theoretically in Rcfs. [35l. [3o|. 

The phase diagram, as a function of the interaction pa- 
rameters g c ( s ) , is contained in Fig. [7] Among the possible 
phases, we find one corresponding to a chiral fixed point 
similar to the case of spinless electrons. In this phase, the 
flow of current is sensitive to the flux through the ring, 
and we find that the charge and spin degrees of freedom 
must circulate with the same chirality. 

We have also found that the inclusion of the spin de- 
gree of freedom allows for the existence of a stable fixed 
point where current flows only between two wires, while 
the third remains uncoupled. Such fixed point was al- 
ways unstable in the case of spinless electrons. Thus, 
in more realistic models that include the electron spin 
into account, it may be possible that controlling small 
anisotropies in Y-junctions of quantum wires may lead 
to sensitive current switches. 
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APPENDIX A: BOUNDARY CONDITIONS AND 
CONDUCTANCE 

We will show in this section how the conductance ten- 
sor is extracted from the boundary conditions. Within 
the linear response theory, we obtain the Kubo formula 
for the conductance tensor of multiple wires introduced 
in Eq.lEUD as in 
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G 3 k,c(s) = lim t r / rfa; 



dTe lUT {T r J jAs) (y,T)J kAs) (x 1 Q)), 



(Al) 



where the currents J c ( s )( x i T ) = — *v2Sr0 c (s) (x, r). The currents can be separated into two chiral currents, Jf~ c ^ = 
\/296j,c(s) = (^x ~ * | 9 r )0 JiC ( s )/v / 2 and Jf c ( s ) = \/2d8j,c(s) = (^x + i9 T )6j tC ^/y/2. In terms of these chiral currents, 
Jj,c(s) = Jf* c (s) ~ Jfcis) anc ^ *he Kubo formula Eg. (|A1[) becomes 



Gjk,c( s ) = un \ I F 



dx [ (T T J* (s) (y,T)Jg c(s) (x,0)) + (T r J^ (A2) 
- (^ c(s) (y,r)j£ c(s) (x,0)) - {T T jf As) {y, t)J r c{s) (x,0))]. 

I 



We will use the rotation matrix 7?. corresponding to the 
boundary conditions in DEBC method to evaluate the 
correlation functions. 



A convenient trick to respect boundary conditions is to 
analytically continue the right mover currents to a; < 
and identify 



1. Boundary conditions and correlation functions 

We first consider the correlation functions of the chiral 
currents in an infinite quantum wire. The off-diagonal 
components {Jf c ( s )Jf c { s )) van ish an d the diagonal com- 
ponents are given by 

(J R (s) (y,r)J R {s) (x,0)) = -2d 2 (9 c{s) (z,z)e c{s) (0))(A3a) 
(Jc( s) (y,T)J^ s) (x,0)) = -2d 2 (9 c{s) (z,z)0 c{s) (Q))(A3b) 

where z = ir + (y — x) . Since the ^-correlation function 



we obtain 



(0 C(S )M)0 C(S )(O)) = -f inM 2 



(J*s)(y,r)J% 8) (x,0)) 



(Jc (s) (y,T)J^ s) (x,o)) = 



z* 

9c(s 



(A4) 

(A5a) 
(A5b) 



Let us recall how the boundary conditions can be writ- 
ten in terms of the rotation matrix 1Z 



$ R {x)=1l T $ L {x)\ x = Q , 
where 6 R a.nd6 L arc defined as 



(A6) 



( <<s) \ 



( C W \ 



and 



(A7) 



for N quantum wires. Because # c ( s ) = y/g c ^ (<f>c(s) 
R \ 



) , the boundary conditions can be translated to 



Jl (s) (0) = ^ s) (0). 



(A8) 



(A9) 



for x > 0. With this identification, the chiral current 
correlation functions between different wires can be eval- 
uated in terms of the matrix elements 



(^c(s)(y,T)J« {s) (x,0)) 



9c(s) x 
9c(s) 



(AlOa) 



(JUM t Ka4 x ^ = -W 5 ^ ( A10b ) 



for the diagonal correlation functions, (RR) and (LL), 
while 



(J*c( S )(y,T)Jj; cis) (x,o)) 



(AlOc) 



(ir+(x + y)) 2 
(ir - (x + y)) 2 
for the off-diagonal correlation functions, (RL) and (LR). 

2. Conductance tensor 



Now, we can insert the correlation functions Eq. (|A10I) 
into the Kubo formula Eq. (|A2[) to evaluate the conduc- 
tance tensor. With the aid of the integral formula, 



dre luJT 



1 



(ir + u)' 



= 2iruj@(u)e- wu , (All) 



where 8 (it) is the Heaviside step function, the Kubo for- 
mula reads 
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G 



ij,c(s) 



hL 



dx [Sij (&{x -y) + Q(y - x)) - H ji Q{x + y) - TZ l3 Q(-x - y)}. 



(A12) 



r 



The integration of the combined first and second 0- 
functions gives a constant L. In addition, since both 
x, y > 0, the integration over the third 0-function gives a 
constant L, while that over the fourth one vanishes. The 
conductance tensor evaluated from the Kubo formula is 
therefore given by 

,:2 



G 



ij,c(s) 



(A13) 



As an example of this generic formula for junctions of 
multiple quantum wires, let us compute the conductance 
for a junction of two wires. First, the N BC is governed by 



the rotation matrix 1Z 



N,ij 



5^. Hence, all elements of 



the conductance Gn vanish, indicating a total decoupled 
junction. Second, the D BC is governed by the rotation 
matrix 



n D = 



1 

1 



(A14) 



Inserting this into Eq. (|A13[ ). one obtains 



Then the modular transformation of the summation in 
Eq. (|B1[) can be achieved by using the Poisson summation 
formula, which replaces the summation by a integration 
with the periodic 6- function, 5s^\(x) 

Z = -Fw / ddx ^A( f ) e_if |S ' 2 ■ (B3) 

■n{q) J 

The S- function can be further written as the sum of A', 
reciprocal lattice of A 



= V7T\ E exP^M* ■ x)} , (B4) 



where Vq (A) is the volume of the unit cell. The partition 
function becomes 



(B5) 



and 



G D -2g c[s) e -( ^ 1 



h = -l2=2g—{V 1 -V 2 ), 
h 



(A15) 



(A16) 



from the definition 7j = GijVj. This is, as one would 
expect, the conductance for perfect transmission in a ID 
quantum wire. 



cast using a standard Gaussian integral. 

To proceed, we shall use the following identity 



f 1 d 

J d d x exp[— — XiAijXj + b ■ x] 



(2ir) d / 2 
Vdct A 



exp[W(b)}, 



(B6) 



APPENDIX B: MULTI-DIMENSIONAL 
MODULAR TRANSFORMATION 

Since modular transformations arc useful in calcu- 
lations using BCFT, we will in this appendix define 
and provide the general formulation of the the multi- 
dimensional modular transformation. 

Generically, a d-dimensional partition function Zq is 
proportional to 



Z = — V q-^ 2 = — V e~T 13,2 



(Bl) 



where q = e~~ and A indicates the lattice points. In or- 
der to check the finite size spectrum with a given bound- 
ary condition, we have to rewrite the partition function 
in terms of q = e~~ , i.e. perform a modular transforma- 
tion. First, the modular transformation of the Dedekind 
^-function reads 



(B2) 



where Aij is a d x d matrix and W(b) = 
\ Ylij=i &i(-^ _1 )»j^3- I n our casc Eq. (|B5|) . one can iden- 
tify that = ^Sij and b = i2ttu' and obtain 



(detA)-, = ( T )- = (^ 



(B7) 



W(b) = - 



(2tt) 2 (3 7T/3_ 



2 2irV 



I 



uT- (B8) 



Inserting these results into Eq. (|B5jl . the partition func- 
tion becomes 



Z = 



V {A)r)(q) d V 2tt/' 



u'eA' 



( 2 ) # v ^ 



(B9) 



where the identity in Eq. (|F32|) is used for second equality. 
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